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We discuss with a rather critical eye the current situation of black hole (BH) solutions in f{R) 
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meaning, existence or lack thereof of a Birkhoff’s theorem in this kind of modified gravity. We focus 
then on the analysis and quest of non-trivial (i.e. hairy) asymptotically flat (AF) BH solutions in 
static and spherically symmetric (SSS) spacetimes in vacuum having the property that the Ricci 
scalar does not vanish identically in the domain of outer communication. To do so, we provide 
and enforce the regularity conditions at the horizon in order to prevent the presence of singular 
solutions there. Specihcally, we consider several classes of f{R) models like those proposed recently 
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the available no-hair theorems, and in the cases where the theorems apply, the absence of hair is 
demonstrated analytically. In the cases where the theorems do not apply, we resort to a numerical 
analysis due to the complexity of the non-linear differential equations. Within that aim, a code to 
solve the equations numerically was built and tested using well know exact solutions. In a future 
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I. INTRODUCTION 


Modified f{R) gravity has become one of the most popular mechanisms to generate a late accelerated expansion 
in the universe without the need of introducing new fields [1]. It was also one of the first consistent models for early 
inflation [2]. During the past fifteen years or so, several specific f{R) models have been thoroughly analyzed in many 
scenarios, but only a few of them can survive the classical tests (e.g. solar system, binary pulsar) while predicting 
the correct accelerating expansion, and in general, a successful cosmological model, both at the background and at 
the perturbative level. Therefore, it is still unclear until what extent this kind of alternative theories of gravity can 
recover all the successes of general relativity (GR) while making new testable predictions. 

As concern black hole solutions, the situation of f{R) gravity can, in some sense, differ from GR, and in other sense 
be almost the same. The last statement is related with the content of Section [ill about the existence of the same 
kind of vacuum black-hole solutions found in GR, while the former concerns the existence of hairy solutions, or lack 
thereof, that we analyze in all the rest of the paper. 

Perhaps the first and simplest theorem concerning BH solutions in GR was the Birkhoff’s theorem (BT). Roughly 
speaking, this theorem establishes that in vacuum all spherically symmetric (SS) spacetimes are also static, and those 
that are asymptotically flat (AF) are represented by a one-parameter family of solutions, namely, the ubiquitous 
Schwarzschild solution, where the parameter is interpreted as the (ADM) mass M of the spacetime (see [5] for a 
discussion). Remarkably, when including an electric field, the SS solution can be extended as to include the charge 
Q of the BH, it is the well known two-parameter Reissner-Nordstrom (RN) solution. In the AF case, both the 
Schwarzschild and the RN solution are i? = 0 solutions of the Einstein’s field equations (also termed Ricci flat 
solutions). When including a cosmological constant A, the Schwarzschild and RN black holes become a two and 
three-parameter family respectively, and the solutions are asymptotically de Sitter (ADS) or asymptotically anti de 
Sitter (AADS), depending if A > 0 or A < 0, respectively. 

In the decade of 1960’s, motivated by the discovery of the Kerr solution, several theorems (notably, the uniqueness 
theorems) were established for stationary and axisymmetric spacetimes both in vacuum and with an electromagnetic 
field (see mu] for details and reviews). One of the main consequences of those theorems is that in the AF case the 
solutions of the Einstein field equations under such symmetries are characterized only by three parameters, the mass 
M, the charge Q, and the angular momentum J of the BH. These solutions are known as the Kerr-Newman family, 
which extends the Schwarzschild and RN black holes to more general spacetimes: stationary and axisymmetric. Due 
to the apparent simplicity of such solutions, Wheeler established the so called no hair conjecture, a statement that 
“doomed” all possible stationary AFBH solutions from having any other parameters than those three (M, Q and J). 
This conjecture has been reinforced thereafter by the elaboration of several no-hair theorems (NHT’s) that forbid the 
existence of BH solutions with more parameters associated with other kinds of matter fields (see [SH] for a review). 
Among such theorems one can mention those that include several kinds of scalar fields. Eventually, this conjecture 
proved to be “false”, for instance, within the Einstein-Yang-Mills system [5], and Einstein-scalar-field system with 
“exotic” potentials that can be negative liiin]. Or when including rotation like in the Einstein-boson-field system m- 
Nonetheless, since most, if not all of such hairy solutions are unstable the community (or at least part of it) consider 
those solutions as weak counterexamples to the Wheeler’s no-hair conjecture. Therefore, it has been tantalizing to 
extend the conjecture in the following more precise, although still informal, statement: the only stable stationary 
AFBH’s are within the Kerr-Newman family [5]. 

The proposal of alternative theories of gravity as a possible solution to the dark-matter and dark-energy problems 
and to other theoretical problems (e.g. inflation, gravity renormalization) has motivated people to generalize several 
of the theorems and conjectures mentioned above, and which pertain to GR, to the realm of other modified-gravity 
proposals. While it is out of the scope of the present paper to review all such attempts, we shall simply focus on 
f{R) metric gravity. Unless otherwise stated, by f{R) gravity we mean a theory that departs from the GR function 
f{R) = R-2A. 

As concerns this kind of theory, a large amount of analysis has been devoted to establish an analogue of the BT 
for the SS situation However, the reality is that no rigorous BT exists today in f{R) gravity, as far as we 

are aware. In fact, if one such theorem were proved, certainly should be restricted to some specific /(i?) models. 
Moreover, the theorem should establish at least four things, upon fixing the boundary conditions (i.e. regularity and 
asymptotic conditions): I) Staticity: the only spherically symmetric solutions in vacuum [i.e. without any matter field 
associated with the standard model of particle physics or any other field that is not associated with the Lagrangian 
f{R)] are necessarily static (i.e. the existence of a static Killing field should be proved from the spherically symmetric 
assumptions); 2) Existence: the existence of an exact static spherically symmetric (SSS) solution in vacuum; 3) 


^ In the sense that a perturbation can lead to an eventual loss of the hair. 
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Uniqueness: the SSS solution found in point ‘2)’ is the only solution in vacuum (or prove otherwise); 4) The conditions 
under which the solution in point ‘2)’ matches or not the exterior solution of an SSS extended body. 

So far, in vacuum only a few BH exact solutions exist in f{R) gravity, and those that are genuine AF, ADS or AADS 
correspond simply to the same kind of solutions found in GR, where R = Rc = const everywhere in the spacetime 
(with Rc = 0, Rc > 0 or Rc < 0, respectively). It is unclear if other solutions exist with the same kind of asymptotics 
but with a varying R in the domain of outer communication of the BH. We shall elaborate more about this point 
below to be more precise. Furthermore, in the presence of matter (i.e. a star-like object), it is possible to find SSS 
solutions where R can vary in space [HHII], however, those solutions are not exact, but given only numerically, and 
it is unclear if the exterior part (i.e. the vacuum part) of those solutions is the same solution found when matter is 
totally absent in the spacetime, if it exists at all, as it happens in GR where the exterior solution of extended objects 
under such symmetries is always given by the vacuum Schwarzschild solution 

Now, despite the absence of such BT, some NHT’s have been proved in this kind of theories. In order to do so, people 
have resorted to the equivalence between a certain class of f{R) models (notably, those where /i? > 0 and /at? > 0, 
where the subindex indicates differentiation) with scalar-tensor theories (STT). The point is that one performs a 
conformal transformation from the original Jordan frame to the so-called Einstein frame where the conformal metric 
appears to be coupled minimally to gravity and a new scalar-field 4>{x) emerges, where x ■= /iti which is also 
coupled minimally to the conformal metric but endowed with an “exotic” potential ^{4>)- Thus, the available NHT’s 
constructed for the Einstein-scalar-field system in GR can be applied for these theories as well (see Section]^, notably 
in vacuum, and when the spacetime is AF and the potential satisfies the condition {(f) > 0 [T^ I20p . 

We stress that the applicability of such NHT’s is possible because the non-minimal coupling between the scalar 
field (f>{x) and the matter fields that usually appears under the Einstein frame obviously vanishes in the absence of the 
matter. The only caveat of this method is that the potential (0) is not given a priori but is the result of the specific 
f{R) model considered ab initio, and thus, ^{(f) can be negative or even not well defined (i.e. it can be multivalued), 
which in turn can jeopardize the use of the NHT’s. Consequently, the existing NHT’s in f{R) gravity can reduce the 
kind of AFSSSBH solutions that are available in some specific models, but do not rule out completely the absence 
of geometric hair. In this context, by (geometric) hairy solutions within f{R) gravity we mean AFSSSBH solutions 
where the Ricci scalar is not trivial (i.e. constant), but rather a function that interpolates non-trivially between the 
horizon and spatial infinity. 

Thus, when the condition {(f) > 0 fails and the NHT’s are not applicable one can resort to a numerical analysis 
for evidence about the existence of such hair or its absence thereof. At this respect it is important to stress that 
regularity conditions have to be imposed at the inner boundary, namely, at the BH horizon r^ in order to prevent 
the presence of singularities there. In Section IV and Appendix]^ we obtain such regularity conditions and then in 
Sections |V A| and |VB| we present analytical and numerical evidence, respectively, showing that hairy solutions are 
absent in several specific f{R) models proposed as dark-energy alternatives in cosmology. In particular, the models 
considered in Section [VB| are precisely those for which the NHT’s cannot be applied as the corresponding potential 
^ {4>) can be negative or is not even well defined. On the other hand, when such hairy solutions are absent, one may 
still find the trivial solution R — Ri = const for which the field equations reduce to the Einstein field equations with 
an effective cosmological constant Agff = Ri/4, and an effective gravitational constant Ges = Go//(i?i) where Ri is 
a solution of an algebraic equation involving f{R) and fn. This includes the case where Ri = 0. Therefore, in such 
circumstances, all the best known BH solutions found in GR exist also in f{R) gravity simply by replacing the usual 
cosmological constant A by Aofj, and the Newton’s gravitational constant Gq by Geff. In view of this we shall argue 
in Section [h] that such solutions are so trivial (i.e. trivial in the context of f{R) gravity) that almost nothing new 
arise from them. 

Finally, we mention that some “non-trivial” exact SSSBH solutions have been reported in the literature as a result 
of very ad hoc f{R) models [121 [21H1S1- Notwithstanding such solutions cannot be considered as hairy solutions 
because they have unusual asymptotics, and therefore, the corresponding “hairless” solution R = const (including 
the Ricci-flat solution) does not even exist with the same kind of asymptotics. We shall discuss one such solution in 


^ In GR and in the stationary and axisymmetric case, the exterior solution of rotating extended bodies does not match exactly the Kerr 
solution. On the other hand, in scalar-tensor theories of gravity (STT) without a potential, and due to the spontaneous scalarization 
phenomenon eg. there exist star-like SSS solutions where the exterior solution is not given by the Schwarzschild metric, since the scalar 
field there is not zero. It turns out, however, that in the absence of ordinary matter, the exterior solution within the same STT is 
given only by the Schwarzschild solution. These two scenarios show that the exterior solutions in presence of ordinary matter are not 
necessarily the same as in complete absence of it, and so this illustrates the difficulty of establishing a generalization of a BT when the 
hypothesis change. 

^ A more recent proof of the same result m adopts a weaker convexity assumption > 0 for the potential. This proof is similar to 

Bekenstein’s m which assumes only stationarity as opposed to staticity. 
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Section IIII Al 


The article is organized as follows: In Section |IT] we discuss in a general setting the conditions for the existence 
of several trivial BH solutions. In Section m we focus on SSS spacetimes and provide the corresponding differential 
equations to find BH solutions. We also discuss some exact solutions that will be used later to test a numerical code 
constructed to solve the equations. The boundary conditions appropriate to solve these equations with the presence 
of a BH are given in Section [m] in form of regularity conditions at the horizon. No-hair theorems and the properties 
of f{R) gravity formulated in the Einstein frame are analyzed in Section]^ In that section we also provide strong 
numerical evidence about the absence of hair for several f{R) models when the NHT’s do not apply. Our conclusions 
and final remarks are presented in Section |VI| Several appendices at the end of the article complement the ideas of 
the main sections. 


II. f{R) THEORY OF GRAVITY 

The general action for a f{R) theory of gravity is given by 


Ildabt'fp] — J 2^ ''Z ~9 ^ 2 :-f/matt [f/ahj ■0] j (1) 

where n = SttGo (we use units where c = 1) and f{R) is a sufficiently smooth (i.e. C^) but otherwise an a priori 
arbitrary function of the Ricci scalar R. The first term corresponds to the modihed gravity action, while the second 
is the usual action for the matter, where represents schematically the matter fields. 

The field equation arising from the action Q under the metric approach is 

fuRab - -^fgab - (VaVb - gab^^) Ir = uTab , (2) 


where fa stands for df/dR (we shall use similar notation for higher derivatives), □ = g‘^'^'^c^d is the covariant 
D’Alambertian and Tab is the energy-momentum tensor of matter resulting from the variation of the matter action 
in Q. It is straightforward, although a non-trivial result, to show that the conservation equation W^Tab = 0 holds 
also in this case (see Appendix [A| for a proof). In turn, this latter leads to the geodesic equation for free-fall particles 
= 0. Therefore, the weak-equivalence principle (for point test particles) is also incorporated in this theory as 
well. Actually f{R) metric gravity preserves all the axioms of GR but the one that assumes that the held equations 
for the metric gab must be of second order. Clearly the only case where this happens is for f{R) = R — 2A, which 
leads to GR plus a cosmological constant (hereafter GRA) 

Now, taking the trace of Eq. ([^ yields 


UR = 


1 

3/fli? - 


KT-3fnRR{^Rf + 2f-Rfn 


where T := T“. When using ([^ in ([^ and after some elementary manipulations we obtain m 


Gab = 


fn - 


fRR^a^bR + fRRR{^aR)i^bR) - ^ {R/r + f + 2 kT') + kT,, 


(3) 

(4) 


Equations Q and ^ are the basic equations that we have used systematically in the past to tackle several problems 
in cosmology and astrophysics [ninMiiss], and that we plan to use in this article as well. 

Now, apart from the GRA theory for which /jj = 1, fuji = 0, and R = 4A — kT, for more general models, one 
imposes the conditions fu > 0, for a positive GeS, and /rr > 0, for stability [52]. However, in this paper we shall 
sometimes relax these two assumptions in order to explore its consequences for the sake of finding BH solutions. 

In vacuum, that is when Tab = 0, or more generally, in the presence of matter fields where T = 0, like in electro¬ 
magnetism or Yang-Mills theory, Eq. ® admits in principle the trivial exact solution R = Rc = const where Rc is 
a solution of the algebraic equation [2f{R) — R/rIZ/rr, = 0. In particular, if Zrr^Rc) yf 0 and 0 < fRR{Rc) < oo. 


4 
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In Section III A 
A result that is 


we extend our units so that Go = 1 as well. 

also a corollary (when applied to a four dimensional spacetime) of a theorem known as Lovelock’s theorem [23]. 
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as often happens in potentially viable f{R) models, then Rc is an algebraic solution of 2f(R) — RJr = 0 that we call 
Ri. For such kind of solutions the field Eq. Q reduces to 

Gab 7 (h) 


where 


Aoff = 


GeS = 


Ri 


Go 


fniRi) ’ 


( 6 ) 

(7) 


we assume 0 < fniRi) < oo. Moreover, the condition 0 < /rrIRi) < oo is an stability condition that is usually 
imposed in order to avoid exponentially growing modes when perturbing around the value R = Ri. 

On the other hand, if there exists other value R — Rc = const, that we call R2 (in order to avoid confusion with 
Ri) such that blows up at i? = i? 2 , but R2 ^ Ri then R2 can be also a possible trivial solution of Eq. 

However, in such a case one must be extremely cautious as in Eq. may appear products of the sort 00 x 0 ^r 
00 / 00 ). Even if such a product is finite, namely zero, still the interpretation of such solution would be problematic 
as Jrr would be singular everywhere. At what extent such “singularity” is physical, is something that one should 
clarify. In Sec. [Ill A[ we shall be dealing with an SSS solution where the Ricci scalar is not constant but R ^ R2 
as r —>■ 00 . Hence, in that example Jr^oo, /rr 00 and Jrrr —>■ 00 as r —>■ 00 . However, these pathologies, as 
peculiar as they may be, does not concern us too much in this article, since we will be mainly interested in situations 
where they are absent. 

Now, in the particular and simplest scenario where R = Ri is one of the trivial solutions of Eq. ([^, we see that the 
field equation (§ corresponds to GR with the usual “bare” cosmological constant A and the bare Newton’s constant 
Gq replaced by Aes and Geff, respectively. Therefore, in that occurrence all the solutions that exist in GR exist in 
f{R) gravity as well when taking into account the above replacements. In particular, the AFBH solutions with R = 0 
that exist in GR, like the Kerr-Newman family and its SSS limit, exist also in f{R) gravity if i?i = 0, i.e., if /(O) = 0. 
On the other hand, the De Sitter or Anti-De Sitter BH solutions associated with the Kerr-Newman family with a 
cosmological constant [301131] exist also in f{R) gravity if i?i ^ 0. 

It is note that BH solutions R — Ri = const, which from the point of view of f{R) gravity are trivial, have been 
systematically reported in the literature as something new or special (e.g. see Refs. uni Eg). However, as we just 
showed, the existence of such solutions stems from the fact that i?i exists in various f{R) models, like the ones we 
consider in Section V A In turn, the existence of such a trivial solution is just a standard demand for f{R) theories 


to produce a late acceleration expansion: the cosmological constant Aeff emerge while the universe evolves towards 
the solution i? —>■ Ri as the universe expands and matter dilutes. So nothing exceptional, astonishing or radically 
different from the already known BH solutions in GR are to be expected in this kind of trivial solutions in f{R) 
gravity. Furthermore, and as a consequence of these remarks, the Bekenstein-Hawking entropy defined for such BH’s, 
has Geff instead of Go in the formula. That is, the modified entropy is S' = A/(4Geff) = AfR{Ri)/{4Go) [33|, where 
A is the area of the BH event horizon, instead of just S = A/ (4Go). 


Let us consider now the “odd” scenario where Ri is such that Jr^Ri) = 0, and Ri = 0. Then Eq. ^ admits i? = 0 
as one possible trivial solution. Nonetheless, the possible solutions for the metric that satisfy Eq. Q degenerate 
and a whole spectrum of solutions can emerge, besides the usual AF vacuum solutions mentioned above. Such 
solutions are all the possible solutions of the Einstein equation Gab = iGTab compatible with a null Ricci scalar R= 0. 
More specifically, if one considers an /(i?) model in vacuum such that 2/ — RJr = 0 at i? = 0 (implying in turn 
/(O) = 0) but assuming /^^(O) ^ 0, then Eq. ([^ is solved trivially. Moreover the field Eq. Q in vacuum reduce 
to /i?(0)Gab = 0 which is trivially satisfied for /h( 0) = 0 even if Gab 7^ 0. Thus, in that instance the model admits 
unexpectedly all the possible solutions associated with the non-vacuum Einstein equation Gab = i^Tab compatible with 
G'jj = —R = —T = 0. Now, given such scenario one has to deal with a problematic interpretation of the new global 
quantities that appear in the BH solutions as integration constants since in reality we are dealing with an f{R) model 
in vacuum. In non-vacuum GR the global quantities (other than M and J) are ascribed to properties associated with 
the matter described by Tab (like the electric charge Q or the nodes of the Yang-Mills field). Therefore the same 
quantities would appear in this degenerate scenario, but clearly they cannot have the same interpretation in vacuum. 
In section [V A| within the framework of Model 1, we shall encounter one such example which is associated with the 
model f{R) = kR^. In this model R = 0 is a trivial solution and clearly /^(O) = 0, thus all possible solutions of 
Einstein equations with a traceless energy-momentum tensor are in principle allowed, for instance, AFSSS solutions. 
In fact, one such solution is exactly the analogue of the RN, where the roll of the charge is played by a new quantity 
that appears as an integration constant. This interesting situation was also analyzed recently in Ref. |34j and remarked 
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before in Ref. [T3]. We shall elaborate more about this in Sections 
uniqueness and the BT. 

Having clarified the fact that f{R) gravity naturally admits trivial BH solutions where R = const everywhere, 
notably the ones with i? = 0, and which correspond to the same solutions found in GR modulo a trivial redefinition 
of constants, the main goal of this paper is an effort to find AFSSS solutions where R{r) is a non-trivial solution of 
the field equations, where r is some radial coordinate. This means that if such solution exists, R{r) should interpolate 
in a non-trivial manner between the event (Killing) horizon of the BH and the asymptotic region. If exist, this is 
what we might call a hairy solution. In order to find such a solution we enforce suitable regularity conditions at the 
horizon. These conditions provide a very specific form of derivatives of several variables at the horizon, notably, the 
first derivative of the Ricci scalar. These regularity conditions are extremely important as they prevent the presence 
of BH solutions that are pathological at the horizon. In particular, they prevent a singularity in the scalar-degree 
of freedom at the horizon, a singularity that would be otherwise considered as physical as opposed to a coordinate 
singularity. In Section |IV| and in Appendix we provide such conditions in its full form and compare them with 
similar regularity conditions reported first in [3S], but amended in [3S]. Examples of such pathologies (singularities) 
are common in STT when using conformal methods to generate SSS exact solutions without enforcing the regularity 
condition of the scalar field at the BH horizon m- Thus, the fact that the field equations may not be satisfied 
at the horizon in those examples cast serious doubts about its relevance as genuine counterexamples to the no-hair 
conjecture. 


HI A and V A notably, concerning the issue of 


III. STATIC AND SPHERICALLY SYMMETRIC VACUUM SOLUTION IN f{R) GRAVITY 


The problem of describing a SSS space-time in vacuum (Tab = 0) within the framework of f(R) gravity reduces to 
solving the field equations Q and (|^ for the following SSS metric: 


ds^ 



2M{r) 

r 


+ (^1 - 


dr'^ -I- -|- sin^ Odip^) , 


( 8 ) 


where the mass function Mir) provides the ADM mass in the asymptotic region provided the spacetime is AF The 
function 5{r) indicates the extent to which the equality is satisfied or infringed by the components of the 

Einstein tensor or equivalently, by the corresponding components of the effective energy-momentum tensor given by 
the r.h.s of Eq. Q taking Tab = 0 there [3H]. So, if d = const, in particular zero, it means that the equality holds 
exactly everywhere in the spacetime. This situation includes some of the the best known SSS spacetimes. The field 
equations for the metric Eq. ([^ have been obtained in a rather convenient form in m with a slightly different but 
equivalent parametrization. We present their final form based on those equations but without the matter terms: 


R” 

M' 




r ) 


1 


(2/ - R/fl)r 


L r-2M 
M ^ 1 

T + 2{2fR + rR'fRR) 


— ^IrrrR'^ 
^IrM 


21 rM' -M) 2 

— 0 - 


(r — 2M)r 


R' 


+ Rpi^fR + /) 
r 3 


rR'fRR\-^ 
3 


Jr 


AM , / 

-Tpi^RfR - f) - Ir + 2ri? JrA 

^ r \ 


1 - 


2M 


r / 


1 


2(2/fl + rR'fRR)r 


2r^ 

— (2/-R/fi) + 


rR'fRR rr^ 


fR L 3 

4(r - + | 

fR 


-{2RfR- f)-2fR 

} 


(9) 


( 10 ) 


( 11 ) 


In order to find BH solutions, Eqs. have to be solved from the BH horizon at r = to the asymptotic 

region which in the AF or Anti-De Sitter cases corresponds to spatial infinity. On the other hand, in a De Sitter 


® Certain spacetimes that are not AF posses a well defined “ADM” mass. In those cases, the Minkowski background is then replaced by a 
suitable background spacetime with respect to which the mass indicates deviations as one approaches the asymptotic region (the mass 
is then part of the “monopole” term in the gu component In those spacetimes the equivalent of the ADM mass corresponds typically 
to a suitable renormalization of M(r oo) (see footnotelsjl. 
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background a cosmological horizon can be reached at some > r^. In this paper we shall focus only in AF solutions, 
therefore we deal only with the BH horizon. 

Suitable boundary conditions at are to be imposed. Typically one assumes regularity conditions that enforce the 
verification of the field equations there (see Section IV). For instance, the value M{rh) = rhl‘2, enforces the existence 
of the Killing horizon, where rh can have, in principle, any arbitrary non-negative value. The condition 5(rh) is also 
arbitrary as it simply determines the value of S{r) in the asymptotic region which can be redefined by a suitable 
change of the t — coordinate. In other words, the field equations are left invariant if one performs the transformation 
(5(r) —)■ (5(r) -I- const, a property that holds due to the existence of a static Killing vector field. So without loss of 
generality we can take S{rh) = 0. As concerns R[rh), hereafter Rh, this value is usually fixed so as to obtain the 
desired asymptotic value for R. Typically, but not necessarily (see below), a value R(rao) = 0 (where Too stands for 
r —>• oo) gives rise to AF spacetimes 

Due to the complexity of the field equations @-([TT]) one usually resorts to a numerical solution, in which case, 
the value Rh is fixed by a shooting method (cf. Section [V B ). This value, which provides the adequate asymptotic 
behavior for R{r), is also related with the the convergence of M{r) to the ADM mass. Equation (10) can be written 
as in GR in the form 


M' = Anr^peff 


( 12 ) 


where pefj can be readoff from Eq. (10). In this kind of coordinates the ADM mass Madm for AF spacetimes is 
given by M{rao), also called the Komar mass [5]. In order for M to converge to Madm one requires peff ^ 
asymptotically with e > 1. That is, we require that the effective energy-density falls off faster than 1/r^. Otherwise 
M(r) can diverge asymptotically as M{r) ~ ln(r) if e = 1 or M(r) ^ r^~’^ if e < 1 ®. As a consequence, if the 
spacetime containing a BH possesses a non-trivial Ricci scalar R{r), this must behave as T^ff ~ —peff, namely 
R ~ with e > 1 in the AF scenario. Therefore, the shooting method is implemented such that R(r) behaves 

asymptotically in the previous manner using Rh as a control parameter. If no such behavior exists for all Rh, this 
implies that AF spacetimes with a non-trivial Ricci scalar do not exist. Thus, scalar-hair would be absent and the 
only possible AF solutions are R — 0. 

There exist spacetimes with asymptotic behaviors different from AF or ADS/AADS with interesting properties. 
One of them correspond to a spacetime that is AF except for a deficit solid angle 0 < A < 1. This kind of spacetimes 
are typically generated by topological defects, like strings and global monopoles m- The deficit angle is related with 
the “symmetry breaking” scale 77 of the Mexican hat potential by A = [4T] . The deficit angle produces typically 

a repulsion instead of an attraction of test bodies in the gravitational field generated by these defects. In this kind 
of spacetimes the mass function M{r) ~ rA/2 asymptotically. It is somehow remarkable that such spacetimes can 
have a well defined mass |42j which, in the case of spherical symmetry, can be computed, in practice, by a suitable 
renormalization of M(roo). This situation corresponds to e = 0. It is then tantalizing to define the “ADM” mass in 
this case simply by Madm A = -^ren(^oo), where M^^{r) = M{r) — ^rl2. In the next section we shall see that in 
fact such a mass is not exactly M^^(roo) but rather proportional to it where the constant of proportionality depends 
on the deficit angle A. Furthermore, R ^ 2Afr'^ asymptotically. Hence, we have here an explicit example where the 
Ricci scalar is not trivial and it also vanishes asymptotically, and yet the spacetime is not AF. Below we present an 
explicit f{R) model for which an exact SSS of this sort can be found, except that the deficit angle is not produced 
by a topological defect, but by the underlying modified-gravity model itself. 


We finish this section by giving explicitly the kind of trivial solution R = Ri = const that we alluded in Section |TT] 
For SSS spacetimes we have 


ds^ = - l- 


2Mo Aefir^ 


dt^ 


dr^ 




' 1 2Mq Apff 

< r 3 


(13) 


^ In Anti-De Sitter backgrounds one demands R{roo) = Ri with Ri < 0. For De Sitter backgrounds enforcing the boundary conditions 
can be more involved since a cosmological horizon at (r^ ^ '^h ^ ‘^) present, and in this case the value Rh is to be fixed so as to 
recover the regularity conditions at as well. 

® In the case of Anti-De Sitter backgrounds the corresponding “ADM” mass is given by Madm = Mren{f'oo) where Mren{f') = M(r)—Ar^/6 
(A < 0). A rigorous definition for the mass in this kind of backgrounds has been given by several authors [Slllo]. This mass is like the 
Komar mass m but specialized to the SSS scenario. For instance, when the mass function is given by M{r) = m + Ar^/6 where m is 
a constant (like in Kottler—Schwarzschild-De Sitter metrics), m coincides with the mass found from the formal definitions |39II4Q| . In 
the context of f{R) gravity peff — -^efT/(^) must behave as with e > 1 for Mren('roo) to converge. 
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with 


M{r) 

= , 

6 

(14) 

Rir) 

= i?l = 4Aeff , 

(15) 

6{r) 

= 0 , 

(16) 

2/(i?i) 

= RifRiRi) ■ 

(17) 


This is the Kottler-Schwarzschild-De Sitter solution, where Mg is a constant of integration of Eq. (101 which is 
identified with the “ADM” mass (cf. Ref. [Snilin])- In the cases where i?i = 0, the solution reduces to the Schwarzschild 
solution. It is a straightforward exercise to check that Eqs. (|l4|)- 0 solve exactly Eqs. regardless of the f{R) 

model, provided the model satisfies Eq. ( |T7| ) and fniRi) ^ 0, ^ 0, fRR{Ri) < oo and fRRR{Ri) < oo. As 

we show in the next section, when the f{R) model is given explicitly Ri writes in terms of its fundamental parameters 
[cf. Eq. (@ ]. 


A. Exact solutions 


Let us consider the model 


/(R) = 2aVR - a = 2a^ 



(18) 


where a > 0 is a parameter with units [distance]”^, and a is another parameter of the model which is related to an 
effective cosmological constant as we show below. This f{R) model and a variant of it was considered in the past by 
several authors [551 [55]. In the SSS scenario the metric 


ds^ 



O 

ar 

IT 


2Q 

J.2 


dt^ + 





(19) 


with the mass function, Ricci scalar, and 5{r) given, respectively, by 


M(r) 

R{r) 

S(r) 


r ar^ 
1 

^ I 

0 , 


Q 

2r ’ 


( 20 ) 

( 21 ) 

( 22 ) 


solve Eqs. (§”([11]) exactly, as one can verify by straightforward substitutions. Here, Q is an integration constant. 
Taking into account electromagnetic and Yang-Mills fields, this solution was extended in Ref. [55]. When a = 0 this 
solution was part of a more general class of solutions associated with the model f{R) = kR^ [T31|53]. However, for 
n = 1/2 the asymptotic behavior of those solutions was not analyzed by those authors as we do here. 

The coordinates are defined such that, — oo <t<oo, O<0<7r, 0<(/)<27r with < r il a < 0 and rh < r < r'/ 
if a > 0, where and corresponds to the location of the event and cosmological horizons, respectively, that we 
analyze below. 

The metric ( |I^ possesses a deficit angle A = 1/2, Madma = 0, “charge” Q and a cosmological constant Aoo := 
i?(roo)/4 = a/^Xsee Appendix [pj for more details). In the current case Madma = (1 ~ [H] where 

= M{r) — Aoor^/6 — rA/2. The divergent terms (linear and cubic in r) appear in this renormalization 
of mass since, as we remarked before, the spacetime has a deficit angle (associated with the linear term) and also a 
cosmological constant (associated with the cubic term). Using (20) we conclude M^^{r) = —Q/ (2r), thus, Madma = 
-Wren(^oo) = 0. Examples of spacetimes with a deficit angle and with zero-mass BH are not new [15] . 

The metric with a deficit angle given by Eq. is a solution with a cosmological constant, Aoo = a/4. For 
instance, taking a > 0 we can introduce = 6/a for convenience. The location of the black hole horizon depends on 
the value of Q. There are three possibilities: a) If Q > 0, Q := q^, the event horizon of the black hole is located at 

Th = \ \j2P- -I- 2P^fl + b) If Q < 0 (Q := —(X) and 1 — 8g^/€^ > 0, there are two horizons which are located 










’’’h — \ _j_ — %(p- jP', and rt = \_ 2PyJl — Sq^ jP. In particular, for an extreme black hole with 

q 2 _ f 2 horizon is given by Vh = |£|/\/2; c) If Q = 0i then = \i\. 

The AADS spacetime with a deficit angle is a solution with a negative cosmological constant A^o = a/4 (with 
a < 0) and -^aads ~ —6/a. Notice that for Q = 0 an horizon does not exist like in the usual Anti-De Sitter (i.e. 
Anti-De Sitter spacetime without the deficit angle). However, for Q = —q^ < 0 an horizon exists and is located at 

When a = 0, the spacetime turns to be AF except for a deficit angle. In this case only for Q = —q^ < 0 there is 
an horizon at = \/‘2q^. Finally when a = 0 = Q, the spacetime is simply the Minkowski spacetime with a deficit 
angle. All other cases correspond to naked singularities. 

In absence of naked singularities a straightforward calculation of other scalars, like RcdR‘^‘^ and RabcdR‘^'’‘^‘^, show 
that the only physical singularity appear at r = 0. Even if a = 0 = Q, such scalars are RcdR‘^‘^ = l/2r'^ and 
RabcdR‘^^‘^‘^ = ^/r'^ while R = 1/r^, and thus, the singularity at r = 0 is entirely due to the deficit angle. Since 
the coordinates used so far do not cover the entire manifold one can look for analytic extensions using Kruskal-like 
coordinates. These extensions and the construction of Penrose diagrams are out of the scope of the current paper and 
will be reported elsewhere. 

Now, for a < 0, there is no cosmological horizon, and thus, one can analyze the solution as r —>■ oo. It is 
then interesting to note that the asymptotic value of the Ricci scalar R{rao) = a is not an algebraic solution of 
2f{R) — Rfn = 0 which is Ri = 4a/3 but rather a pole of /rh = —a[R — a]“^/^/2. The trivial solution i?i = 4a/3 
provides a cosmological constant Aes = i?i/4 = a/3 which is different from the actual one Aoo = a/4. Now, even 
if fn = a[R — Irr = —a[R — a]“^/^/2 and /rrr = 3a[i? — a]“®/^/4 blow up at i? = a, the quantities that 

appear in Eq. ® behave well asymptotically (i.e. they are finite) as i? —>■ a: 


fRRRi^R)^ 

Irr 

2/ - RfR 
Irr 


fRRRg^^R'^ g''^R'^ ^ 

fRR R- a ^ 2r2 

2a 

-2(R-a)(3i?-4a) ~ ^ 


(23) 

(24) 


Moreover the quantities that appear in the r.h.s of Eq. Q behave also well asymptotically and give rise to the 
cosmological constant Aoo = a/4 as we show next. Take for instance the r — r component of Eq. Q in vacuum. For 
our purposes it is more convenient to take the mixed components. Then 


r 9 


(iRR^r^rR , fRRRR'^\ I f„ , / ^ 


Now let us analyze the asymptotic behavior of each of the terms at the r.h.s of Eq. (25): 

9rr ' 


12 


g 

fRR^ rR 
Ir 


ar^ 

2 

ar 
12 ’ 

V V /? 

- / , =- Vr^rR 

2{R-a) 2 ^ 

( 6 


/ 6 

ar 

2 lr4 

12 


grrR 


IrrrR' 

/fl 

f 

R+^ 

Jr 


2 \r 
3R'^ 


6 ag, 

.4 


6j.2 

3 

4(i? — a)^ 


— R 2(R — a) — a A 


6 


9 ’ 


Therefore, to leading order we obtain 




a 

6 


a a 
12 ~ 6 


a 

'4 


(25) 

(26) 
(27) 


(28) 

(29) 

(30) 

(31) 


9 





















so the r.h.s of Eq. Q is well behaved asymptotically, and it is just a constant which we can precisely identify with 
effective cosmological constant Aoo = a/4, with a < 0. Notice that this constant emerged not only from the last two 
terms at the r.h.s of Eq. (251 but also from the contribution of the first two, which one would naively think that they 


do not contribute as i? —?► a asymptotically. However, as mentioned in Sec. |TTj a closer look shows that one actually 
has in those two terms something like oo x 0 asymptotically. This is why it was necessary to perform the correct 
asymptotic analysis which leads then to contribution —a/12 due to the first two terms of the r.h.s. of Eq. (251. Of 
course one can perform the same asymptotic analysis in the full set of equations -([TT]) to find that all of them 
behave well and consistently as i? —)■ a, since from both the l.h.s and the r.h.s one obtains exactly the same behavior. 
This is otherwise expected as we have explicitly the exact solution from which one can compute M' and R' and R" 
to confirm that nothing diverges as i? —)■ a. The definition of this cosmological constant Aoo = a/4 is consistent with 
the canonical form that the metric coefficients gtt and take in ( |19[ ) in these coordinates. For instance, in terms of 
Aoo they read gu = l/ffr-r = -(1 - A - where A = 1/2. 

Hence, we conclude that when fRu has a pole precisely at the Anti-De Sitter point, the cosmological constant Aoo 
does not arise simply from the last term of Eq. (25), like in the analysis performed in Sec. [TT| where Aeff = i?i/4. That 
analysis was valid provided that as i? —>■ i?i the following two necessary conditions were satisfied: Ri = 2/(i?i)///j(i?i) 
and < oo, which as emphasized above, is not the actual case for this exact solution. Finally, we mention that 

for a = 0, which corresponds to a null cosmological constant, the quantity (2/ — RJr)/Jrr = —6i?^ that appears 
in Eq. (1^ vanishes as the solution approaches the asymptotic value i? = 0, even if Jrr —>■ oo. Moreover, the term 
fRRRR/fRR = —3i?'^/(2i?) also vanishes asymptotically since i?' ^ 1/r^ and R ~ l/r^. Therefore, a posteriori one 
can understand why Eq. ([^ is well behaved asymptotically. 

It is important to stress that in the previous works [H121 US] the above physical and geometric interpretation of 
the metric (19) was completely absent and therefore, its meaning was rather unclear. In the nonvacuum case, some 


but not all of the aspects discussed above for this exact solution were elucidated [ 2 ^ . 

So far we have mainly discuss the case a < 0. Regarding a > 0, the Ricci scalar will tend to a asymptotically but 
will never reach this value since well before r —>■ oo the cosmological horizon is reached by the solution. Therefore, 
i?'(r£) ^ 0 and R"{r’^) ^ 0. That is, the possible solution i? —>• a with i?' = 0 = R" is never reached asymptotically 
due to the presence of the cosmological horizon. 


Now, as concerns the trivial solution R = Ri = const, the model (18) admits the solution i?i = 4a/3, which solves 
i?i = 2/(i?i)//i^(i?i). Notice that = —a[a/3]“^/^/2. Therefore the SSS solution is given by the metric (13) 

where 


Aoff = 


Ri 


a 

3 


(32) 


In this case the cosmological constant is given by Agff, instead of A^o = a/4. 

We emphasize that the solution © is not approached asymptotically by the solution ( |19[), s ince as we mentioned 
before, the latter has a deficit angle A = 1/2, while (13) has A = 0. In summary, the model (18) admits the two exact 
solutions (19) and (13), one where the Ricci scalar is constant everywhere and one where it varies with the radial 


coordinate r. However, it is important to remark that both solutions are not two different solutions with the same 
boundary conditions, but two different solutions with different boundary conditions. Even if we put Mq = 0, Q = 0, 
one solution still has a deficit angle while the other does not. Moreover, both solutions have different cosmological 
constants Aeg and Aoo, and so the spacetimes are not even the same asymptotically when a 7 ^ 0. On the other hand, 
one could try to make coincide the inner boundaries (i.e. the event horizons) of both solutions artificially as well as 
the value Rh by fixing Q and Mq, but still both spacetimes would have different global quantities (mass and charge) 
besides the respective deficit angles (A = 0 and A = 1 / 2 ) and the respective cosmological constants (Agff and Aqo). 
Namely, one solution would have zero “ADM” mass and non zero charge, while the other solution would have non-zero 
“ADM” mass and zero charge. 


Two final remarks are in order. The first one concerns the deficit angle solution (19). It is possible to show that 
in fact one can obtain the /(i?) model (18) by using a kind of “reconstruction method”. This consists by imposing 
a solution in the form of Eqs. (19) with (20)-([2^. That is with g^r = —1/gtt and with M and R written in terms 


of finite powers of r. The condition grr = —^/gtt leads then to S{r) = const, i.e., S'{r) = 0, which in turn provide a 
differential equation for f{R) in terms of r, R{r) and M(r). However, given such a power expansion, one can invert and 
write r = r{R), to obtain a differential equation for f{R) and R solely, which when solved provides (18). This kind of 
“tricks”, which are also common to find exact solutions in GR with exotic sources [H] , have been applied in modified 
theories of gravity in the past [25l[26]. Physically, one usually proceeds in the opposite way we just described. That 
is, one constructs or proposes an explicit f(R) model in order to fit some observations, for instance, in cosmology, and 
then one asks if such a model admits or not an exact solution in one or other scenario. At this regard, and as a second 
remark related with our comment at the end of the previous section, we mention that almost all physically viable 
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f{R) models admit the trivial solutions R = Ri = const given by the Schwarzschild-De Sitter/Anti-De Sitter solution 
Eqs. (13|-(16), where the parameter a is replaced by a different constant according to the different f(R) models. For 
instance, in Section [V A| we will see more specific examples. These are the kind of trivial solutions that we alluded 
in Section |n] that several authors have systematically reported as something “special” in f(R) gravity, while we see 
that they are nothing more than the usual solution found in GR with A —>• Aos = i?i/4 and Go —>■ GeS = Go/fii(Ri)- 
In particular the mass parameter Mq can be redefined as Mg = GesM where A4 can be taken as the fiducial mass 
associated with the spacetime. Moreover, the RN or the Kerr-Newman BH solutions with or without Agff can also be 
found in f(R) gravity when considering matter with a traceless energy-momentum tensor as these corresponds also 
to the same trivial solutions R = Ri = const, including the AF ones when Ri = 0. 

We remark that there are more exact solutions for other choices of n in the model f{R) = kK^ OEl], not only for 
n = 1/2, which corresponds to the model we have just analyzed taking a = 0. Nevertheless, most of those solutions 
are still trivial or have exotic asymptotics. In Section V A (cf. Model 1) we discuss the case n = 2, which was not 
covered in Ref. Ham], as their equations become singular precisely for n = 2. 


Finding physically interesting exact BH solutions different from the trivial ones proves to be difficult when the f{R) 
model is complicated, like the physically viable models that have passed many cosmological and Solar-System tests 
(e.g. Models 4-5 of Section VA|. In particular, if one is interested in genuine AF or ADS/AADS type of spacetimes. 
In such an instance, one has then to appeal to a numerical analysis. At this regard, we stress that we have used the two 
exact solutions presented in this section as a testbed for a FORTRAN code developed to solve numerically Fqs. (§- 
(11) for more complicated and “realistic” /(i?) models, submitted to suitable boundary (regularity) conditions that 

when an horizon is present). These regularity conditions are presented 


represent the presence of a black hole 
next. 


i.e. 


IV. REGULARITY CONDITIONS 

In order to obtain the regularity conditions at the horizon r = r^j, whether is the event (inner) horizon or the 
cosmological (outer) horizon, we expand the variables as follows: 

F{r) = F{rH) + {r - rH)F'{rh) + l:{r- rH)F”(rH) + “ rH)^F'"{rH) + 0{r - rn)'^ (33) 

2 o 

where F{r) stands for M(r), i?(r), 5(r). When replacing these expansions in Eq. and demanding that the 
derivatives of these variables are finite at the horizon one obtains after long but straightforward algebra the following 
regularity condition: 


R'\ 


2r{RfR-2f]fn 


r2(2R/fl - /) - 6fR 


Jrr 


(34) 


where as stressed, all the quantities in this equation are to be evaluated at the horizon r = r^, which in principle is 
any positive value, whether rn is the event (inner) horizon (denoted by rh) or the (outermost) cosmological horizon 
(r/). In this paper we will only be interested in finding AF solutions, and thus, the spacetime will contain only the 
event (inner) horizon. 

In turn, Fq. (101 provides: 


M'l 


r-2 (^2RfR - 


IVr 


(35) 


Hereafter the quantities evaluated at the horizon will be written with a subindex 


^ If one takes = 0, the regularity conditions correspond to a spacetime with a regular origin. In this case, R'\r=o = 0, M\r=o = 0, and 
M'\r=o = 0. The so called solitons are localized field configurations that are globally regular^ including the origin. In GR there exist 
examples of spacetimes with matter that allow for solitons and SSSAFBH with hair, like in the Einstein—Yang—Mills system EllIS] and 
in the Einstein-scalar-field system with an asymmetric scalar-field potential 
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It is to be mentioned that regularity conditions of this kind were proposed first in |35j . and then rectihed in [36j . 
Those authors did not use exactly the same differential equations (101 and 0 that we used here, which were a 
consequence of Eqs. and® . Instead, they departed from the much more involved field equation (H). Nevertheless 
our regularity conditions (IM I and (35) are equivalent to those of [35] . 


In order to illustrate the consistency of these two regularity conditions, we take three f{R) models and their exact 
solutions. First, f{R) = R — 2A, which corresponds to GR plus a cosmological constant. The SSS solution is like 


in (131 but with A instead of A^ff. In this case Eq. (34), when multiplied by gives an identity 0 = 0, since the 


numerator at the r.h.s R/r — 2/ = i? — 4A = 0. On the other hand, since i? = 4A = const then R' = 0. Therefore, 
R' Jrr = 0. On the other hand, the r.h.s of Eq. (35) yields M'^ = Ar|j/2, which corresponds precisely to the De 
Sitter/Anti-De Sitter value obtained from Eq. (20). 


The second case corresponds to the trivial constant solution discussed in previous sections: R{r) = i?i = const, 
i?'(r) = 0, where i?i given by {R/r - 2f)R^ = 0 and such that fRR{Ri) ^ 0, < oo and fRRR{Ri) < oo. 

Then Eq. (34) gives = 0, which is clearly compatible with the trivial solution, whereas the r.h.s of Eq. (|3^ yields 


M'jj = Riff^S, where i?i = 4Aeff. 

A less trivial test to our regularity conditions is provided by the exact solution given by Eqs. (20)-(22). For that 
solution it is easy to verify that both sides of Eq. (34) give = —2/r^ while Eq. (35) yields 1/4 + ar'^18 + 

Q/{4rjj) in both sides. 


Now, Eqs. ( [^- ( 11) do not provide the regularity conditions for i?" and <5' at the horizon. We require to differentiate 
Eqs. j9) and ( |11[ ) with respect to r one more time. When doing so and when replacing the expansions in the form 
Eq. (1^) it is possible to obtain the values R'^ and 5'fj. This is a lengthy but otherwise straightforward calculation. 
In the Appendix |B] we provide the explicit expressions. 


V. NO HAIR THEOREMS AND SCALAR-TENSOR APPROACH TO f{R) THEORIES. 


As remarked in the Introduction, there are very well known NHT’s for the Einstein-scalar held system, hereafter 
Einstein-^ system (i.e. Einstein-Hilbert gravity minimally coupled to a real scalar held </>) when the potential associated 
with a scalar held 4> verihes the non-negativity condition {(f)) > 0. On the other hand, it is also well known that 
provided Jrr > 0 one can map f{R) into a Scalar-Tensor Theory (STT) by dehning the scalar y = /r. This is the 
so called Jordan-frame representation of the f{R) theory. Under this framework and in the presence of matter, the 
scalar held y turns out to be coupled non minimally to the curvature but minimally to matter. A further conformal 
transformation allows one to write the theory in the so called Einstein frame, where a new scalar held (f> with a new 
potential ‘^ (</>), is now coupled minimally to the curvature but non-minimally to the matter sector. Nonetheless, in 
the absence of matter, and still under the Einstein frame, the theory look exactly as an Einstein-c/ system, and thus, 
one can turn to the NHT’s and see whether they apply or not for the f{R) model at hand. Since we are interested 
precisely in hnding SSS black holes in vacuum, we can thus exploit this equivalence to see if hair is absent or if it 
can exist. Thus, given a specihc f{R) model, we carry out the following protocol: 1) We check if the model allows 
a “trivial” solution R = 0 (leading to Schwarzschild solution). If the model does, then it means that the model is a 
priori able of allowing AF hairy solutions (i.e. solutions where R would be a non-trivial function of r that interpolates 
from the event horizon, with value Rh, to spatial infinity where R = 0). Many f{R) models used as geometric dark- 
energy admit also the trivial solution R — 0. Thus, in this paper we focus only in such kind of models and analyze 
if geometric hair can also exist. In a future work, we shall analyze the existence of hair in ADS or AADS black-hole 
solutions; 2) We write the model as a STT in the Einstein frame. If the potential (</>) is not negative, then the 
NHT’s apply and we conclude that AFSSS hairy solutions cannot exist in such model. In fact, if ^{(f>) turns to be 
strictly positive, the Schwarzschild solution does not even exist, only the De Sitter type of solutions exist at best; 3) 
If the potential has negative branches then we proceed to analyze if hairy solutions exist by solving numerically the 

IV and the Appendix [B| 


field equations of Sec. HI under the regularity conditions provided in Secs. 


In order to fix the ideas, let us now briefly review the different scalar-field transformations required to analyze 
properly an f{R) model under the STT approach. 

Let us consider for simplicity the following gravitational action without matter 




^grav — / d ~g 


fQmR-Q) + f{Q) 


(36) 
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where at this point Q is an auxiliary scalar-field that depends on R and fq = df /dQ Now, if fqqiQ) ^ 0 (in 
particular fqqiQ) > 0) the variation of this action with respect to the metric give rise to field equations equivalent 
to the original theory 0 in vacuum provided Q = R. Thus the action ( |M| is equivalent to 0 if we enforce Q = R 
(see Appendix [C] for a further discussion). 

Furthermore, one can introduce the scalar field y := fqiQ) and write a dynamically equivalent action as 


djF =Jd'^x X - X^V(x)J, 


(37) 


where V(x) is defined as follows: 


V(x) = 


2kx^ L 


Q(x)x- f(Q(x)) 


(38) 


By dynamically equivalent^ we mean that we obtain field equations that are completely equivalent to the original 
action (see Appendix [C|). To do so, we can treat the scalar-field x metric-independent, and R as independent of y, 
and thus take the action as a functional of both gab and y. Thus, variation of this action with respect to the metric 
leads to field equations equivalent to Eq. ([^, while variation with respect to y simply leads to R = Q. We see then 
that the action ( [37| ) is equivalent to the action of a Brans-Dicke like theory, with a Brans-Dicke parameter wbd = 0 
(implying that the kinetic term associated with the gradients of y is absent) and with a potential 17(y) = 2Ky^E(y) 
(see Appendixfor the details). Clearly a STT of this sort (i.e. one with wbd = 0) but without a potential would 
be incompatible with the bound wbd ^ 4 x 10^ [46], which is required for the theory to pass the Solar System tests, 
and thus, would be automatically ruled out. However, the presence of this potential makes possible for certain f{R) 
models to pass those tests even if wbd = 0. This depends if the potential allows for the emergence of the chameleon 
mechanism (cf. Ref. |47l|48]), but not all the potentials have this property. 

While we will consider some /(i?) models that seem to pass such tests thanks to the chameleon mechanism, our 
main purpose in this paper is to analyze the issue about the existence or absence of hairy BH solutions in such models 
rather to test their observational viability. Thus, there are models that we use only for that purpose and which may 
fail the observational tests. 

It is said that, the actions (36) and (37) are written in the so called Jordan frame (JF), where by frame is to 


be understood as a set of physical variables. As we remarked, in reality we have only introduced a new scalar-field 
variable without changing the metric in any way, thus this formulation of the theory is completely equivalent to the 
original one (at least in the sectors where fan > 0), and thus, the field equations are also equivalent. 

As mentioned before, the energy-momentum tensor of matter is conserved in the JF (see the Appendix [A|, and 
therefore it is with respect to this frame that point-test particles follow geodesics. 

Furthermore one can define a new scalar field and a conformal metric as follows: 




9ab — X 9ab — ^ 


= py ^ 


9ab 


In terms of these new variables the gravitational action (37) takes the form 


IeF =ld*X - ^i4>) 


(39) 

(40) 

(41) 


where all the quantities with a tilde are defined with respect to the metric gab- This corresponds precisely to the 
Einstein-Hilbert action coupled minimally to a scalar field cj). Notice that, unlike the JE action, a kinetic term appears 
due to the conformal transformation. The action (41) is written in what is known as the Einstein frame (EE), where 
the potential is defined as ^{(f) = 17(y[(/']), which will be given explicitly when the model /(i?) be provided 

(see Section [V A ). If we include the matter action then the scalar field (/> will be coupled non minimally to the matter 
fields. In this paper we are only interested in the vacuum case so the field equations obtained from the action (41) 


Do not confuse the field Q with the parameter appearing in the metric ||l9|. 


tion § 


we have omitted in the integral the term y ^ V“Va.'?5>, which can be converted into a surface term, that we assume 


In the action 
to vanish. 

See Ref. |55| for a thorough discussion of STT in the EF with the presence of matter. 
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are simply those of the Einstein-(/> system: 


Gab 

rpCf) 

^ab 

Ucj) 


kt: 


ab ’ 


(Va<?i>)(Vf,(/>) - gab 




d4> 


(42) 

(43) 

(44) 


As we emphasized previously, based on this equivalence and in view of our interest in finding SSS and AF non¬ 
trivial black holes in f{R) gravity, we have to take into account the NHT’s [El[20] which are valid when {(/)) > 0 
The theorems roughly establish that whenever the condition ^ {(/)) > 0 holds, given an AFSSS spacetime containing a 
hlack-hole (with a regular horizon) within the Einstein-(j) system, the only possible solution is the hairless Schwarzschild 
solution. Here by hairless we mean that the scalar held (j){r) = ^O; he., the scalar held is constant everywhere in 
the domain of outer communication of the BH and it is such that = 0, in order to prevent the presence of a 

cosmological constant which would spoil the AF condition. 

The NHT’s can be avoided if the potential has negative branches, notably at the horizon illE]- So in our case, 
given an f{R) model, we have only to check if the corresponding potential satishes or not the condition > 0. 

In the affirmative case, we conclude that SSS and AF hairy black holes are absent. Nevertheless, when this condition 
fails, one usually need to resort to a numerical treatment in order to analyze if a black hole can support scalar hair 
or not. 

Before concluding this section, a final remark is in order. There is an important relationship between the critical 
points of the potential notably the extrema, and those of the “potential” y{R) defined such that [T7](see also 

Appendix [C| 

w) = % = liy-RfR), ( 45 ) 

In Section |h] we denoted the extrema of F(R) by i?i. It is straightforward to verify 


^ ( 2f-RfR \ ^ 

d(l) \d(l)J\dx) k/| V2k/ V 3 / \2k) IndR 


(46) 


Therefore, provided fn^Ri) ^ 0 and fn^Ri) < oo (i.e. the conditions for a well defined map to the Einstein frame), 
we see that the extrema of correspond precisely to Ri. However, care must be taken when /^{Ri) = 0 or 

fniRi) = 00 , as it may happen in some models that we will encounter in the next section. When this happens, the 
conformal transformation ( |40| ) becomes singular or ill defined as x = Ir- 

Finally, for /(i?) models where fun is not strictly positive, notably, where fun can vanish at some R = Rw, called 
weak singularity (cf. Model 5 in Section V A) it will be useful to introduce the “potential” W (R) defined via 


= -IR = 


2/ - RJa rdm 


S/rr 


Irr 


(47) 


The finite or divergent behavior of Wr{Rw) also provides insight about 1^(7?^,). For instance, if Wr^Rw) is finite, it 
means that yR{Rw) vanishes like fRR{Rw)- Furthermore, Wr{R) can suppl y furth er information about the possible 
trivial solutions R = const. As we discussed in Section [TT] and also in Section HI A if at i? = i ?2 where //j_R(i? 2 ) = oo 
and #fi(i? 2 ) vanishes, R = R 2 can be one trivial solution of Eq. (|^ in vacuum or more generally, when the matter 
has a traceless energy-momentum tensor, a solution that can be different from R = Ri if Yr{R 2 ) 7 ^ 0. 


A. f{R) models and the NHT’s 


In this section we focus on some f{R) models that also admit the trivial solution R = 0, and check whether or not 
they satisfy the condition ^ (cj)) > 0 for which the NHT’s apply. The results are summarized in Table |l] at the end of 
this section. In order to obtain the EF potential ^ {4>) from V (x) let us recall the relationships 


X 

^^(0) 


vix[4>\) , 


(48) 

(49) 


The theorems actually account for multiple scalar fields. 
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the first one is obtained from Eq. (39), while the second one is a definition. 


( \ n 

^ j , where n, A„ and are positive parameters of the theory. fixes the scale for each 
n and A„ can be chosen to be proportional to i?„. This model has been thoroughly analyzed in the past in several 
scenarios (see [55] and references therein). For this model, 


nA„ / i? \"-i 


y=fj^ = A") 

A J n \ D I 

^ 71 . ^ -^n ' 


(50) 


If we focus in the domain R S [0, oo), then for n > 1, y € [0, oo). In fact the value y = 0 corresponds to a degenerate 
situation Geff —>■ oo that we discuss below. For 0<n<l, y— ^ooasi?— and vice versa. So, y € (0, oo) for 
0 < n < 1. 


By inverting the relationship (50) and using Eq. (38) followed by the use of Eqs. (48) and (49), the two potentials 
read 


v{x) = 




2Kn VnA„ 


X 


2-n 

n-1 




defined for y > 0, and 




(51) 


(52) 


defined for — oo < (j) < +oo. The condition '^(^) > 0 holds if n > 1, while for n < 1, the potential is ^(</>) < 0, but 
in both cases the potential does not have minima, at least not for a finite </>. Moreover, we do not consider n < 0 
because /a < 0 which can give rise to an effective negative gravitational constant. For n = 1, corresponding to GR, 
V{x) = 0 = '^(0) as one can see directly from Eq. (38). For this model the NHT’s a priori apply since > 0. 


The fact that the potential is strictly positive and have no minima imply that the solution R = 0 cannot even exist 
in the EF. The point is that for 0 < n < 1 the solution i? = 0 corresponds to y —)■ oo {(j> —>■ oo), while for n > 1 the 
same solution corresponds to y = 0 = i? (</> —)■ —oo). We see then that in both cases the mapping to the STT in the 
EF is ill defined precisely at y = 0 where gab = 0, while —>■ 0 as y —?> 00 This problem at y = 0 exacerbates 
for n = 2 that we discuss below. We conclude that for this f{R) model AFSSS simply cannot exist under the EF. In 
the original formulation, the theory also degenerates at i? = 0 for n > 1 since /i?(0) = 0. AFSSS solutions exist but 
they are no unique as we are about to see. 

As we remarked briefly at the end of Section |llj for this class of f{R) models a quite degenerate situation may 
occur. To fix the ideas, let us focus on the case n = 2 in the original formulation, since in the STT approach the 
maps breakdown at i? = 0 as we just mentioned, given that X = Ir = const x R. So in this case the field Eqs. (|^ 
and (§ in vacuum reduce to 


Gab = ^Wa^bR-^R^ 
K I 4 . 

□i? = 0 . 


(53) 

(54) 


where we used /rrr = 0, (R/r + f)/fR = 3R/2 and (2/ - RfR)/fRR = 0 in Eqs. (53) and (54) 


Therefore we see that i? = 0 is a trivial solution of Eq. (54). On the other hand, for such trivial solution Eq. (53) 
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is satisfied for any Gab ^ 0 compatible with G\ = —R = 0. For instance, this can be satisfied for any solution of 
the metric satisfying the Einstein equation Gab = RTab^ with T“ = 0. This degeneracy is somehow remarkable as 
shows that solutions of the field equations in f{R) gravity may not be unique, as illustrated by this simple model. In 
the AFSSS scenario, one BH solution is clearly the Schwarzschild solution, but other solutions are possible [nisi], 
which would be interesting to know to what kind of matter content they correspond in pure GR. In this example, the 
AFSSSBH solutions are unique as concerns the solution R — 0, since at the horizon R = 0 = R', but they are highly 
non unique as concerns the metric. That is, to the trivial solution ii = 0 of Eq. ([^, one can associate any solution 
for the metric satisfying Gab = ^GTab, with T“ = 0, like the Schwarzschild solution, the RN solution, a solution within 
the Einstein-Yang-Mills system, etc. 


This particular model illustrates the care that one has to take when transforming the original variables to the STT counterpart. In the 
next section we shall encounter situations where (</>) is not even well defined as it turns to be multivalued. 

For this f{R) model (2/ — RfR)/fRR = {ttfR + /)//fl = Thus we can expect degenerate solutions (in the sense 

described in the main text) for n > 2 as well. 
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As emphasized in [HIM], this f{R) model provides a specific example showing that a generalization of the Birkhoff’s 
theorem similar to the one elucidated in the Introduction, simply cannot exist in general. It is enlightening to stress 
that this degenerate situation in vacuum is in a way similar, but opposite, to the “degeneracy” that appears in GR 
with matter sources: given f{R) = R, Eq. (§ or Eq. @ reduce to the Einstein field equation, whereas (§ reduce to 
/rrOR = 0, with /rr = 0. This means that this equation is satisfied identically regardless of the value □!?, which in 
general, is not zero because R = —kT. In other words, in GR the metric is constrained to satisfy the Einstein’s field 
equation, but R is not constrained to satisfy any differential equation like Eq. ©■ 


Finally we mention that when R = const = i?i, the model wi th n = 2 also admits the trivial solution Eqs. (13)- 
( |l7l ) Incidentally, for this model the algebraic condition ( (T7| ) is satisfied for any Ri. That is, i?i emerges as an 
integration constant independent of the parameters of the model. Therefore the value Aes = Ri/i depends on the 
assigned value for Ri. In particular, taking i?i = 0 we just recover the usual Schwarzschild solution as mentioned 
above. 

For any other n > 1 the model admits the trivial solution i? = 0, but the solution for the metric is not unique 
either as the degeneracy emerge as well in a similar way to the case n = 2 . 

Model 2: f{R) = R + C 2 Ri{R/Ri)'^, where C 2 is a positive dimensionless constant and Rj is a positive parameter 
that fixes the scale. This model was proposed by Starobinsky as an alternative to explain the early inflationary period 
of the universe [2]. For this model. 


X = /fl = 1 + 2c2i?/R/ 


(55) 


In principle the model is defined for —oo < R < oo. However, if we impose X = /jj > 0, then we require R > 
—Rj/{2c2). In fact if we allow x ^ 0 the transformation to STT is not well defined, and the model degenerates at 
X = 0 in the original variables as GeS —>- 00 . If we focus on solutions with R > 0 then x ^ 1- Furthermore, Jr and 
Irr are both finite at i? = 0 for any C 2 S K. 

Proceeding like in the previous model, this one has associated the following potentials 


V{x) 


Ri 

8c2K 

Ri I 
8c2K 


\ - e~ 


(56) 

(57) 


The potential ^{(p) is defined for —oo < <j) < oo. The region (p > 0 corresponds to y > 1, while —oo < (p < 0 
corresponds to 0 < y < 1 . 

Glearly ‘^{(p) > 0, with a global minimum located at ^ = 0 where the potential vanishes (see Figure]^. The 
model admits the solution R — Q corresponding to y = I and (p = 0. The trivial solution i? = 0 is the only root of 
l^(i?) = R/3. For this model the NHT’s apply, and therefore AFSSSBH hairy solutions are absent. 

Although the model is defined for C 2 > 0 in order to be compatible with inflation, in the context of BH’s one can 
in principle consider C 2 < 0 as a way to evade the NHT’s because then (cp) < 0 (in that instance the condition 
y > 0, implies R < i?//( 2 |c 2 |) and the global minimum becomes a global maximum). It turns out, however, that 
potentials that are negative around a maximum but vanishes there (in this case ( 0 ) = 0 ) sometimes admit non¬ 
trivial solutions that seem hairy and AF^^. We shall discuss a numerical example of this kind later, but suffice is to 
say that the BH solutions that one finds may vanish asymptotically but can have an oscillatory behavior that make 
them not genuinely AF. In order to illustrate this, suppose that asymptotically the metric component grr behaves 
like Qrr ~ 1 -|- 2 C'ir'^ sin(C 27 ’)/r (where cr and Ci ^2 are some constants, and 0 < tr < 1 ), then the mass function 
M{r) ^ sin(C' 2 r) oscillates (it may even diverge if cr ^ 0), and thus, it does not really converge to a finite value 
in the limit r —> oo, a value that one would identify with the ADM mass. Yet g^r —>■ 1 as r —> oo. Thus, for this kind of 
solutions the spacetime is not authentically AF. These arguments can be justified using the following heuristic analysis. 
Let us consider Eq. (44) and neglect the non-flat spacetime contributions from the metric. Moreover, expanding {(p) 
around its maximum (which is equivalent to expand —^{(p) around its minimum) gives ‘^{(p) = —m‘^(p'^/2. With 
these simplifying assumptions it is easy to see that the SSS solution of Eq. (44) is 


(p = (po 


sin(a; -I- Xq) 


(58) 


For n = 2 the potential ^ (0) is constant in the EF, thus the trivial solution 0 = const gives rise to the Schwarzschild-De Sitter solution 
just like in the original variables. 

In more general scenarios such potentials lead to instabilities because perturbations around the maximum leads to growing modes. Thus 
the trivial solution 0 = 0 would be unstable and never settle into a stationary configuration. 
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FIG. 1: (color online). Potential (in units of Ri/Go\ cj) is given in units of associated with the Model 2 with 

C 2 = 1. The potential verifies > 0 for any C 2 > 0 and therefore, the NHT’s apply. The minimum at i)) = 0 is associated 

with the Schwarzschild solution where R{r) = 0. 


where x = mr, and xq, </>o are constants. Now, at leading order when r —>■ oo, we can take the conformal factor 


X = 1 in Eq. (40) and both metrics (the Jordan and the Einstein frame metrics) coincide asymptotically. Thus the 


energy-density contribution = —T^ ’ ^ is given by 


^ j .>2 . </>§ cos[2(a;-h xo)] 

P<P=2^ = - - —2 -- 


(59) 


which is not positive definite. As a consequence, the mass functions is not positive definite either. In fact, the mass 
function behaves asymptotically as M{r) ^ (/)Qsin[2(a; -I- aio)]/m, which oscillates with r and does not converges to 
a definite value (the Komar mass). As we mentioned above, this heuristic analysis confirms the behavior of the full 
numerical solution for potentials of this kind. One such sort of solutions are shown for the Model 4 in Section |VB| 
below. 


Model 3: /(i?) = R — Oii?* In , where a\ is a dimensionless constant, and i?* is a positive parameter 

that fixes the scale. This model was proposed by Miranda el al. [16] as a possible solution for the accelerated 
expansion of the universe while being free of singularities during the cosmic evolution. This model is also viable for 
constructing relativistic extended objects. Notwithstanding, it has problems at the level of the Solar System and 
grow of structure [49j. Although it may not be a realistic model in all the scenarios, it is worth consider it in this 
context due to its simplicity. The model is defined for R > —i?* and the condition /i? > 0 restricts R in the range 
R > i?*(ai — 1). In particular for R^, and ai positive, as we assumed for this model, whenever /^ > 0 is satisfied the 
condition R > —R^ is also satisfied. Furthermore fan = aii?y^(i?/i?* -I- 1)“^ is finite at A = 0, as well as / and Jr. 
For this model 


X = /i? = 1 


CTli?* 

R^ R 


(60) 


The domain i?*(ai — 1) < i? < oo corresponds to 0 < y < 1, and to —oo < (j) <0. The potentials are 


V{x) 


i?* 

2kx^ 


ai In 


ai 

1-X 


+ 1-X-OLI , 


A, 

2k 




f 

ai \ 


1-eVW‘^J 



+ 1 - Qfi 


(61) 

(62) 


The potential '^(</>) is thus defined for —oo < ^ < 0. 

The solution i? = 0 corresponds to y = 1 — ai in this model. However, for ai > 1 one is led to y < 0, which 
by construction is not allowed in the EE frame. Therefore this solution cannot be present in that frame for those 
values of ai. In fact, for such values of ai the potential {4>) has a minimum at some but > 0 

(see Figure]^ According to our protocol, a strictly positive potential cannot allow for the Schwarzschild solution. 


The minimum at 


corresponds to one of the roots Ri 0 of which leads to a Schwarzschild—De Sitter type of solution. 


17 







































FIG. 2: (color online). Left panel: Potential (in units of Rt/Go', <(> is given in units of associated with the 

Model 3 for different values of a\ in the range 1 < ai. For such values is strictly positive. In particular for oi = 1 the 

potential vanishes when (j) —^ —oo. For ai > 1 the potential has a global minimum at <()i,min which leads to a De Sitter type of 
solutions <()(r) = (f> 2 ,inin in the EF and R{r) = const in the JF. Middle panel: Potential ‘i/ (cf)) for five values of ai in the interval 
0 < ai < 1. For 0 < ai < 1, the potential satisfies > 0 and it vanishes at the global minimum 02,min leading to the 

Schwarzschild solution where 0(r) = 02,min in the EF and R{r) = 0 in the JF (see the right panel). Right panel: the function 
dW/dR is depicted. The zeros of this function at i? = 0, which are mapped to 02,min in the EF, lead to the Schwarzschild 
solution. The zeros at i? = — 1 (in units of i?*) in principle lead to a Schwarzschild-Anti-De Sitter solution, but / —>■ oo and 
fn ->• -oo, /rr ->• oo, /rrr ->■ -oo there. 


Therefore, such solution cannot be recovered from the EF approach. As a consequence, the NHT’s also rule out the 
existence of AFSSS hairy BH solutions. In the original variables the trivial solution can be recovered since R — 0 is 
a root of but then the effective gravitational constant Gee becomes negative, as the condition /_r > 0 fails at 

i? = 0. 

Now, for ai = 1, (0) —> 0 if 0 —> —oo, which is the minimum, a situation similar to the Model 1. Finally, for 

0 < «! < 1, the potential satisfies '^(0) > 0 (see the middle panel of Figure]^. In this case the NHT’s also apply, 
and the Schwarzschild solution in the EF is associated with a local minimum 0 = 02,min at which ^(02,min) = 0. 
The trivial solution 0(r) = 02,min is associated with the trivial solution i? = 0. 

The overall conclusion for this model is that hair is in general forbidden, i.e., non-trivial AFSSSBH solutions R{r) 
attempting to interpolate between Rh and R — 0 cannot exist if fa > 0. More specifically: 1) For ai > 1 the EF 
representation precludes the presence of hair in AFSSSBH according to the NHT’s since the solution i? = 0 does not 
even exist as the potential (0) never vanishes (the potential is only defined for /_r > 0). If in the original JF variables 
one permits the possibility of having negative values for //j (which in turn implies negative values for Geff), then the 
trivial solution R = 0 may exist. But then when looking for a hairy solution that interpolates from Rh to i? = 0, the 
solution may cross the value = 0 if JrIr^ > 0, a value which is associated with the singularity GeS = Gq//r oo. 
Such hairy solution would be rather pathological if it exists at all. On the other hand, if we impose the condition 
/r < 0, in order to avoid that singularity, then R turns to be restricted in the range —i?* < R < i?*(ai — 1). In 
particular one would require ai > 1 to allow for the solution i? = 0 to exist. This range for R restricts severely the 
region in which one can look for an optimal shooting Rh (see Section VB|, unless 1 <C ai. All in all, the numerical 
exploration shows that hairy solutions seem to be absent in all these scenarios; 2) For 0 < ai < 1 the NHT’s apply 
straightforwardly, and therefore hairy solutions cannot exist either. Unlike the previous subclass (ai > 1), the solution 
R{r) — 0 corresponding to 0(r) = 02,mm (02,min) = 0] can also be recovered in the EF since for this range of values 
of Q!i the condition fR{0) > 0 is fulfilled. 


We turn now our attention to three f{R) models that have been analyzed recently in cosmology and which are some 
of the most successful ones as concerns the cosmological and the Solar-System tests. However, as we show below, for 
two of these models (i.e. Models 5 and 6) the potential (0) is not even well defined as it turns to be multivalued. 
Therefore, the conclusion about the absence or existence of hair in these two models has to be obtained numerically 
using the original formulation (Sections [TT| and HI) as opposed to the EFSTT approach. 


Model 4: f{R)=R — ReXe{l — e «=) where Ag is a dimensionless constant that usually is taken to be positive for 
a successful phenomenology and Re is a positive parameter that fixes the scale. For instance, Ae > 0 ensures /rr > 0 
and /rr never vanishes at a finite R regardless of the sign of Ag. This exponential model has been analyzed in the 
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past by several authors [iziisg. For this model the scalar field x is 

X = /fl = 1 - . (63) 

The condition //j > 0 holds provided R > i?elnAe, in which case 0 < x < 1, where x 1 as i? —)■ oo. In particular, 
for 0 < Ae < 1 the condition /^ > 0 is satisfied if i? > —i?e|lnAe|. The potential ^{R) has only a global minimum 
at i? = 0, i.e., 1^(0) = 0 [57], and thus, it leads to the Schwarzschild solution (cf. the right panel of Figurej^. For 
Ae > 1 the potential 'f'(R) has a local maximum at i? = 0, and the potential develops in addition a local minimum 
at some R < 0 and a global minimum at some R > 0 |27|. These extrema correspond to trivial solutions R = Ri 
associated with the Schwarzschild-Anti-De Sitter and Schwarzschild-De Sitter solutions respectively. Notice, however, 
that /i?(0) < 0 since the condition R > RelnXe fails at i? = 0. Therefore, in this case the Schwarzschild solution has 
GeS < 0. Since the mapping x —>■'(' is defined only for x > 0, we conclude that for Ae > 1 the condition R > R^lnXe 
implies R > 0 for the EFSTT approach to be well defined and thus, like in the previous model, one cannot recover 
the solution R — 0. As a consequence, we require the original formulation of the theory to analyze if hairy solution 
can exist for these values of Ag. But, again, a hairy solution, if exists, can encounter the singularity at fn = 0 before 
approaching the asymptotic value /^(O) < 0, notably if JrIr^ > 0. 

The potentials are 


V{x) 


Re 
2rx^ L 


(X - 1) ln( 


1-X 


) + Ae + X — 1 



+ Ae + - 1 


(64) 

(65) 


which are valid in the domain 0 < x < 1 —oo < 0 < 0, respectively. The potential > 0 is depicted in 

Figure]^ (left panel) for various values of Ae > 0. For this model the NHT’s apply. As emphasized above, for Ae > 1 
the potential (c^) is strictly positive, in particular its minimum, and therefore the corresponding solution R = 0 
cannot be recovered from the EFSTT approach, but only the De Sitter type of solution R = Ri = const which is 
associated with the minimum of '^ (</>). On the other hand, for 0 < Ae < 1 the potential vanishes at its minimum 
which leads to the solution R = 0 that allows one to recover the Schwarzschild solution. 

If for a moment we dismiss the condition Jrr > 0 and consider Ae < 0 we can evade the NHT’s because now ((/>), 
defined for 0 < (f> < +oo, is never positive (see the middle panel of Figure]^. The new domain is a consequence 
of Eq. (63) which yields x = 1 + |Ae|e~^. Thus for — oo < i? < +oo the scalar field x is defined in the domain 
1 < X < + 00 , which in turns leads to 0 < ^ < +oo. The potentials depicted in Figure (middle panel) suggest 
that hairy solutions might exist. For instance, in the JF variables, such solution R{r) would interpolate between the 
horizon Rh and its value i? = 0 at spatial infinity. It turns out that indeed such solution can be found numerically (see 
the next section) but the spacetime is not authentically asymptotically flat as discussed above within the framework 
of Model 2. The the mass M{rao) never converges to a well defined value as the mass function behaves asymptotically 
M{r) ^ r'^g{r), where 0 < cr < I and g{r) is an oscillating but presumably a bounded function (cf. Figure]^. 


Model 5: f{R) = R + XrRs 




-9 


- 1 


where As is a dimensionless constant, q a dimensionless parameter 


and Rs is a positive parameter that fixes the scale. This model was proposed by Starobinsky m as a mechanism for 
generating the late accelerating expansion while satisfying several local observational tests. We analyzed this and the 
Model 3 in the past in the cosmological setting |5^ and for constructing star-like objects m using the approach of 
Sectionl^ In this paper we take As = I and explore several values of q (see Section VB). For this model the conditions 
fR> 0 ^r a positive Geff) and /rr > 0 do not hold in general. In fact, Jrr vanishes at i? = R^ = ±Rs/y/2q + 1. 
Since /rr appears in the denominator of Eq. ( [^, the vanishing of Jrr was termed by Starobinsky a weak singularity. 
One can appreciate these features from EigureH (right panel) where the “potential” W(R) is depicted. We see that 
\Wr{R)\ = oo at Rf where Jrr vanishes. Thus, the weak singularities at Rf cannot be “cured” by the term 2f — R/r 
in YrIR) because such term does not vanish there, and which otherwise could have lead to a finite Wr. Therefore 
any solution R{r) intending to interpolate between Rh and R — 0 such that Rh > R^ or Rh < will irremediably 
encounter the weak singularities at i?J where we expect a singular behavior in Eq. (pi) . As a consequenc e, ou r search 


for a numerical BH solutions with non-trivial R was limited mostly in the range \R^ < (see Section VB) 


>19 


In the cosmological scenario one usually aims at a De Sitter “point” Ri 7 ^ 0 (as opposed to the Minkowski “point” Ri = 0) in order 
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FIG. 3: (color online). Left-panel: Potential (0) (in units of Re/Go; (j) is given in units of associated with the Model 4 

depicted for three values of Ae > 0. The potential satisfies (0) > 0 for 0 < Ae and therefore the NHT’s apply. In particular for 
1 < Ae the potential is strictly positive and has a global minimum which leads to a De Sitter type of solution R = Ri = const. 
For 0 < Ae < 1 the potential can vanish at the global minimum which is associated with a solution R = 0. The panel displays 
one such example for Ae = 0.5 (cf. right panel). For Ae = 1 the potential vanishes when (j) —>■ —oo. Middle-panel: Potential 
^ {(f)) depicted for three values of Ae < 0. We appreciate that {(j>) < 0 and thus the NHT’s can be evaded. The potential 
vanishes at the global maximum, and therefore the solution R = 0 exists (see the right panel). In principle the potential 
suggests that a hairy BH solution (fi{r) might exist interpolating between some (ph at the horizon and (jimax at infinity where 
'^(<(>max) = 0. However, the numerical analysis shows (cf. Figure]^ that such solutions are not genuinely asymptotically flat 
as the mass function at infinity is not well defined (i.e. it never converges to a hxed value). The local minimum ‘^((('min) < 0 
that appears for each Ae < 0 is associated with the Anti-De Sitter solutions (see right panel). Right panel: the function 
d'y{R)/dR for four values of Ae shows the zeros at which the trivial solutions R — Ri = const exist, including the solution 
R = 0 (corresponding to the places where ^ {(f>) vanishes as depicted in the left and middle panels) and the Anti-De Sitter 
solutions (one zero i?i < 0 for each Ae < 0 is appreciated) corresponding to (j>min where the local minimum ^{(fimin) < 0 as 
shown in the middle panel. 


For this model the potential 'V{R) has several extrema (see the middle panel of Figure]^, in particular, a global 
minimum at i? = 0 with fniO) = 1, fRR{0) = —2Xsq/Rs, which allows one to recover the Schwarzschild solu¬ 
tion. Notice that the global minimum at i? = 0 corresponds to the global maximum of W{R). This is because 
d^W/dR^\R=Q = /dR?\R^Q and d‘^'f'/dR‘^\R=o is positive, whereas /^^(O) < 0, and so d'^W/dR‘^\R=o is neg¬ 

ative. On the other hand Yr{R) at is well behaved there. The other extrema, a local maximum and minimum, 
lead to two Schwarzschild-De Sitter solutions with positive R = const. 

Now, the inversion R = R{x) required to recover the potential V (x) and then the potential ^ ((/>) demands Jrr > 0 
or Jrr < 0. That is, the inversion is possible when x — /fl is a monotonic function of R, which is not the case 
for this model. In principle one could perform the inversion picewise in very specific domains of the model but not 
in all the domain where the model is defined. In view of this drawback the potential {(j)) is not well defined. In 
fact it is multivalued as we are about to see. Its expression cannot be given in closed form but only in parametric 
representation through the equations 


x{R) 

, ^ 2Xsq{R/Rs) 

JR -L p 11+g ’ 

1 + (R/Rs)^ 

(66) 

HR) 


(67) 

^(HR)) 

:= V{x[HR)]) ■ 

(68) 


The form of the potential ^{(f)) is shown in Figure]^ (left panel). Given that ^{(f)) is not single valued it is a priori 
unclear how to establish a method to solve the differential equations in the EFSTT approach and decide unambiguously 
which value of ^ (0) to assign for a given (f). Hence we conclude that one cannot obtain any rigorous result from 


to recover an effective cosmological constant asymptotically (in time), and thus, to mimic the dark energy. In that scenario the actual 
numerical solution R(t) is always positive and larger than thus, the solution never crosses the weak singularity [56]. Something 
similar takes place for the Model 6. 
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FIG. 4: (color online). Left panel: Potential (in units of Rs/Go', <l> is given in units of associated with the 

Starobinsky Model 5 with As = 1 and q = 2. The potential {(p) is multivalued and has negative branches, therefore, the 
net’s cannot be applied. The arrows indicate the trajectory of the parametric plot for increasing values of R {R is in units 
of Rs)- One of the marks (in gray) indicate the value R — 0 corresponding to 0 = 0 at which = 0. The second mark (in 
blue) indicates the starting point of the parametric plot at R = —10. Incidentally, for i? —>■ oo the field 0 —>■ 0 [cf. Eq. ( |66| 
and the potential returns to its starting point. Middle panel: Potential T'(_R) showing explicitly the extrema (maximum or 
minimum) where the trivial solutions R{r) = const exist. The global minimum at 7? = 0 leads to the Schwarzschild solution 
with R{r) = 0. The other extrema are associated with the Schwarzschild-De Sitter solutions R{r) = const. For this particular 
model any non-trivial solution R(r) interpolating between Rh (the value at the BE horizon) and R = 0 (the asymptotic value) 
is confined within the range —1/\/5 < R < 1/y/^ which corresponds to \R\ < 0.447, i.e., values of R close to the global minimum 
at R — 0. Outside this range, the weak singularities can be reached by R (see the right panel). Right panel: the function 
W{R) is depicted showing the places where Wr{R) = dW/dR diverges. These places called weak singularities are located at 
R± = ±l/\/5, where /rr = 0 (these values are denoted generically by R 2 in the main text). At such values Eq. ^ blows up. 


this frame using this potential, let alone trying to implement the NHT’s. But even if we tried to do so, the lower 
branch of the potential does not satisfy the condition ^ {(p) > 0 required by the theorem to prevent the existence of 
hair. In view of this, any strong conclusion about the existence or absence of hair must be obtained from the original 
formulation of the theory that was presented in Section |III[ Furthermore, due to the complexity of the model itself 
and of the differential equations, a numerical analysis is in order. In the next section we provide the numerical results 
that show evidence about the absence of hairy AFSSS black holes in this model. 


I R ) 

Model 6: f{R) = R — i?HS — / ^ " —> where ci and C 2 are two dimensionless constants, and like in previous 

models, i?HS fixes the scale. This model was proposed by Hu and Sawicky [47], and it is perhaps one of the most 
thoroughly studied f{R) models (see Ref. |S5] for a review). In the cosmological context, Ci and C 2 were fixed as 
to obtain adequate cosmological observables, like the actual dark and matter content in the universe. For instance 
taking n = 4, their values are Ci « 1.25 x 10“^, and C 2 « 6.56 x 10“® (56). Notice that the Model 5 with <7=1 and 
this model with n = 2 are essentially the same. Like in the previous model, the conditions fR>0 and /rr > 0 are 
not met in general, therefore, the potential ^ {(p) is multivalued and has negative branches as well. It can be plotted 
using a parametric representation as in the Model 5: 


x{R) 

m) 

^{cpm 


, , nci(i?/RHs)”-^ 

/fl = f - -p - — 

1 + C2(i?/RHS)" 

Vixm]) . 


(69) 

(70) 

(71) 


Figure depicts the potential {(p) (left panel) where one can appreciate the pathological features. In fact, in this 
model a weak singularity /rr = 0 is located precisely at R — 0, i.e., the value that R{r) should reach asymptotically 
in the AF scenario, and it is also the value corresponding to the (hairless) Schwarzschild solution that we should be 
able to recover. Nevertheless, and unlike Model 5, this singularity in Eq. ^ or in Eq. © disappears for some values 
of n because Wr is finite or vanishes at i? = 0 in this model. Namely, Wr vanishes at i? = 0 for 0 < n < 2. Thus, in 
this range of n, the model admits the trivial solution i? = 0. We do not consider the case n = 0 as this model reduces 
to f{R) = R + const, which amounts to GR plus a cosmological constant. For n = 3 we find Wr{0) = —R^g/{18ci), 
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FIG. 5: (color online). Left panel: Potential ‘^((/>) (in units of i?Hs/Go; (f) is given in nnits of associated with the 

Hu-Sawicky Model 6 for n = 4 with the values ci and C 2 as in the main text. The potential {(j>) is mnltivalned and has 

negative branches and the NHT’s cannot be applied. The arrows and marks have the same meaning as in Figure]^ Middle 
panel: Potential fF(_R) showing explicitly the extrema where the trivial solutions R{r) = const may exist. At the local minimum 
R = Q occurs a weak singularity where fan = 0 (see the right panel). Right panel: Potential 'W{R). Two weak singularities 
where Wr(R) = oo are located at i? = 0 and R « 2.34. In particular, the singularity at i? = 0 precludes the search for numerical 

AF hairy solutions with R{r) —>■ 0 as r —>■ oo since the “singularity” is encountered at finite r. 


which does not even vanish. Therefore, this means that R = 0 does not solve Eq. Q, not trivially nor asymptotically. 
For n > 2 and n 3, there is indeed a wieafcsingularity at i? = 0 where \Wii\ = oo (3. the right panel of FigurejS). For 


n < 1 the quantity fniO) blows up, thus we consider only n > 1, notably, for the numerical analysis of Section VB 
In this model the potential iF(i?) has an minimum at i? = 0 for any n > 0, which allows for the Schwarzschi 
solution whenever fjm ^ 0 at R — 0 (i.e. for \n\ < 2). However, when fuR = 0 at i? = 0 and for which |#r( 0)| = 
oo the Schwarzschild solution may not even exist. In those situations, non-trivial solutions where R{r) vanishes 
asymptotically will encounter such singularity (see Section V BI. 

Like in the Model 5, any analysis using the ill-defined potential {(j>) for the Hu-Sawicky model is not robust. We 
then turn to a numerical analysis using the original formulation of the theory. This is presented in the next section. 


B. Numerical analysis and the qnest for hairy solutions 


As we discussed in the previous section, in some circumstances it is possible to formulate the original f{R) model 
as a STT in the EF where the scalar-field turns to be coupled minimally to the EF metric but it is subject to a 
potential '^{</>). If this potential verifies the condition > 0, then the NHT’s apply and, at least in the region 

where Jr > 0 and fun > 0, we can assert that hair (where (j){r) or equivalently R{r) are not trivial solutions) is 
absent, in which case, the only possible AF solutions are at best (j>{r) = const and R{r) = 0. This conclusion follows 
for the Models 1-4 in the sectors where their parameters allow for the R = 0 solution and led to ‘W {(p) > 0. On the 
other hand, we mentioned that Models 2 and 4 can have potentials {p) with negative branches if we allow for the 
parameters C 2 and Ae to be negative. Negative values of such parameters are not usually considered in cosmology, 
but for the sake of finding hairy solutions, we can contemplate them. Because the NHT’s do not apply when ^ (p) 
is negative, notably at the horizon, the problem of hair reopens when this happens. At this regard, several strategies 
are available to solve it: 1) Show an explicitly exact AFSSS black hole solution with hair; 2) Prove analytically the 
absence of it (i.e. extend the NHT’s); 3) Show numerical evidence about one or the other. 

Given that the differential equations presented in Section |HI| are very involved, strategies 1 or 2 might lead to a 
dead end, thus we opted for option three. In particular, this strategy seems even the most adequate as concerns the 
Models 5 and 6, where the potential {p) is not even well defined. 

We proceed to solve numerically Eqs. (§-([n]) subject to the regularity conditions at the horizon provided in Sec.|IV| 
and in the Appendixj^ The only free conditions are the value rt and Rh- The methodology is roughly as follows. One 
starts by fixing the size of the black hole r^, and then looks for R^ so that i? —?► 0 as r —-boo. This “boundary-value” 
problem is solved using a shooting method [53] within a Runge-Kutta algorithm. We have implemented a similar 
methodology for constructing star-like objects in f{R) gravity in the past [17] . Numerical solutions with non-trivial 
hair with asymmetric (non positive definite) potentials have been found previously within the Einstein-(/) system using 
similar techniques [Sj. As we will see in the next section, for certain /(i?) models is not even necessary to perform a 
shooting as the dynamics of R naturally drives i? —)■ 0 asymptotically for a given Rh- 
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Model 

Potential in EF 

Properties 

No-hair 

theorems 

/(R) = A„(£)" 



Q , > 0 , An > 0 , n > 0 

/ 

/(R)=R + C2(^)' 


>0 ,C2 > 0 

'^(0) <0 ,C2 < 0 

/ 

X 

f{R) =R-aiR^ ln(^l-b 


-b 1 - ai 

‘^(</>) >0 ,!<«!, R. > 0 

>0 , 0<ai<l , R. >0 

/ 

/ 

f{R) = R-ReXe{l-e-^) 


- l) 

-1- - 1 

>0 ,1 < Ae, Re > 0 

'^(0) >0 ,0<Ae<l, Re>0 

<0 ,Ae < 0, Re > 0 

/ 

/ 

X 

f {R) = R XsRs ~ ^ 

(ill defined: multivalued) 

- 

- 


(ill defined: multivalued) 

- 

- 


TABLE I: f{R) models and their corresponding scalar-field potentials in the Einstein frame. The last column indicates if 
the no-hair theorems apply (/) or not (X). The potentials that are strictly positive definite > 0, do not even admit 

Schwarzschild BH’s with R{r) = 0. Models 5 and 6 lead to potentials ‘i/ (</>) that are generically multivalued and which have 
no definite sign. In those models the EFSTT approach is not well defined and the applicability of the no-hair theorems is in 
jeopardy, thus, we perform a numerical analysis in the original variables in order to find evidence about the existence or absence 
of geometric hair in AFSSS black holes. 


Now, as we mentioned previously, for the AF solutions to exist, it is not sufficient that i? —)■ 0 as r —>■ -boo. In 
Section [ill A| we analyzed one exact solution where this happens precisely, and yet, the solution is not AF but has a 
deficit angle. In that case the mass function M(r) diverges at least linearly with r. It is then crucial to ensure that 
the mass function converges to a constant value (that we assume to be the Komar or, equivalently, the ADM mass) 
in order to claim for a genuinely AF solution. 

As a matter of fact, we used that exact solution as a testbed for our code. That is, we took the model Eq. (18) as 
mput and recovered numerically the exact solution provided by Eqs. notably for a = 0. Notice that in this 

case i?(r) is not trivial. Eigure § depicts the analytic and the numerical solutions superposed, showing an excellent 
agreement between the two. Typical numerical errors are depicted in Figure 

We also checked that the trivial solutions R = Ri = const that exist in several of the Models 1-5 were recovered 
numerically when starting with R^ = Ri and which lead to the hairless Kottler-Schwarzschild-De Sitter solutions, 
including the plain AF Schwarzschild solution when i?i = 0. 

Additionally we devised other internal tests to verify the consistency of our code. These tests are similar to those 
implemented in our analysis of star-like objects m, and are independent of the fact that exact solutions are available 
or not. 

Let us turn our attention to the specific models that deserved a detailed numerical exploration. 


Model 4: We consider the Model 4 with Ae < 0. In this sector of Ae the potential {(j>) is not positive definite 
and the model may admit hairy solutions because the NHT’s do not apply. Notwithstanding, the only solutions with 
a non-trivial Ricci scalar R{r) that we find numerically are not exactly AF. The Ricci scalar vanishes asymptotically 
in an oscillating fashion as r —^ -boo, but the mass function M{r) does not converge but oscillates as well and grows 
unboundedly as r°' with 0 < cr < 1 (see Figure]^. This behavior is similar to the one provided by the heuristic 
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FIG. 6: (color online). Left-panel: Ricci scalar R{r) (the exact and numerical solutions are superposed) computed using the 
model Eq. (181 for a = 0 (i.e. null cosmological constant). The Ricci scalar is not trivial and vanishes asymptotically, however, 
the spacetime is not exactly asymptotically flat but has a deficit angle [see Eqs. ( |19[ )-(|2^ ]. At the horizon the Ricci scalar 
satisfies the regularity conditions. Middle-panel: the mass function M(r) is not constant but grows linearly with the coordinate 
r due to the deficit angle. Right-panel: metric components —gtt, Qrr and their product —gtt x grr = In the middle and 

right panels the exact and numerical solutions are superposed as well (cf. Figure]^. 



Logio(r/rh) 


FIG. 7: (color online). Relative error between the numerical and the exact solution for R{r) as depicted in Figure]^ (left panel). 
Similar relative errors (not depicted) are found for the mass function M(r) and the metric components. 


analysis within the Model 2, except that here we take into account the full system of equations. Despite such behavior 
the metric components, which depend on M{r), remain bounded in the asymptotic region. This can be partially 
understood by looking to grr = 1 — 2M(r)/r, and realize that for Vh r the non-oscillating part of this metric 
component behaves like So if ct < 1, that part of grr may converge to 1 very slowly, so slowly that one cannot 

even notice it by looking to the numerical outcome. 

This behavior seems to be generic for any Ae < 0 and any Rh- The conclusion is that we do not find any genuinely 
AFSSS black hole solution in this model. 


Finally, let us focus on the Models 5 and 6 that led to pathological potentials in the EFSTT description. 


Model 5: For the Starobinsky model we limit our search for a shooting value Rh first in the region 0 < i? < 
Rs/\/2q + 1 and then in —Rs/y/2q~+T < i?/, < 0 in order to avoid crossing the weak singularities at ±Rs/\/2q~+T 
when R{r) tries to reach the asymptotic value R = 0. We never found a successful shooting parameter leading to 
an authentic asymptotically flat solution. Two examples of this kind of solutions are depicted in Figures and [TO 
Figure]^ shows that for g = 2 the solutions are similar to the exponential Model 4 with Ae = —3 depicted in Figure]8 
Thus, the asymptotic behavior does not correspond to an AF spacetime. 

For g = 4, we find situations where the Ricci scalar decreases monotonically to a constant value without oscillating as 
one can see in the left panel of Figure[^ However, this constant is not related with the trivial solution R = Ri = const 
which is the solution of the algebraic equation 2/(i?i) — i?i/fl(i?i) = 0. In fact, what happens is that M{r) —oo as 
r —>• oo, as we can see from the middle panel of Figure 10 and also M{r)/r —oo, therefore, by looking at Eq. (|^, 
we appreciate that the combination [2/(i?) — i?/fl;]/(l — 2M/r) goes to zero even if 2/(i?) — Rf^ ^ 0. In this way, we 
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FIG. 8: (color online). Numerical solutions for the Model 4 taking Ae = —3. Left-panel: Ricci scalar for three different values of 
Rh- The solutions vanish asymptotically. Middle panel: the mass function M(r) corresponds to the solution of R with Rh = 1 
shown in the left panel. Similar plots for M(r) are found for the other two solutions of R. The non-converging behavior of 
M(r) as r —> -|-oo indicates that the spacetime is not AF. Right panel: metric components for the solution depicted in the 
middle panel. The metric components —gu and Qrr and their product —gtt x grr = are bounded but oscillate as r —^ -|-oo 
corroborating that the resulting spacetime is not AF. 


see that i?' —?► 0 and R" 0, while R const asymptotically, which solves Eq. d^. This behavior of M{r) explains 
why the metric components vanish asymptotically (see the right panel of Figure ^lO| . Before vanishing we see that 
= 1 at logiQ{r/rh) « 0.6 precisely where M{r) = 0. At this value of r, the component gtt = —Thus, we 
conclude that the AF behavior is not recovered either for this and other values q > 0 and different Rh- 

In summary, we find strong numerical evidence that for the Starobinsky Model 5 AFSSSBH with geometric hair do 
no exist. This conclusion is obtained by changing the parameters in several combinations as well as the values Rh- 

Model 6: The numerical analysis of the Hu-Sawicky Model 6 requires more care because fun. vanishes at i? = 0 
for several values of n (|n| > 2). The weak singularity where Wr{R) = (2/ — Rfn)/J rr diverges can be reached 
where Jrr = 0, except if 2/ — RJr vanishes at the same R- In fact for |n| > 2 and n ^ 3 the quantity #r( 0) always 
diverges. In Figure]^ (right panel), one can appreciate this divergence for n = 4. For n = 3 the quantity Wr{0) 
remains finite, but Wr{0) 7 ^ 0, so i? = 0 cannot be a possible asymptotic solution of Eq. ([^. Thus, the case n = 3 
is irrelevant for AF solutions. In consequence, for \ri\ > 2 and n ^ 3, the weak singularity is approached as R tries 
to reach its asymptotic value and at this point Eq. (1^ becomes singular. Therefore, in this range of n we never find 
a well behaved solution regardless of the shooting parameter Rh and R systematically shows a divergent behavior in 
the numerical solutions. 

In the interval 1 < n < 2 the quantity #^(0) vanishes and we do not find any pathologies as i? —>■ 0. The numerical 
solutions are very similar to the oscillating solutions found in the Starobinsky Model 5 and the exponential Model 4. 
Figure shows a prototype of such solutions. As we remarked before, for n = 2 the Hu-Sawicky model is essentially 
the Starobinsky model with q = 1, thus it not surprising to find, at least for these values of n, such an oscillating 
behavior for all of our trial values Rh- Our conclusions seem to be insensitive for several values of the constants Ci 
and C 2 - 

Hence, the numerical evidence indicates that AF solutions with non-trivial hair are absent as well in the Hu-Sawicky 
model. 


VI. CONCLUSIONS 

We have argued that generically f{R) gravity contains trivial solutions R = const- This property allows to mimic 
an effective cosmological constant that can produce the observed accelerated expansion in the universe. In the 
context of black-holes (stationary and axisymmetric, or static and spherically symmetric) the same property allows 
to find the same solutions known in GR (under the same symmetries) with or without a cosmological constant. The 
only difference between both type of solutions (i.e. the solutions found in one or the other theory) is that the two 
fundamental constants involved (the Newton’s gravitational constant and the cosmological constant) are replaced by 
the effective ones, GeS and Agff in f{R) gravity. Therefore, all solutions of this kind reported in the literature within 
the framework of f{R) gravity do not provide any deeper knowledge than the ones we already know in GR. 

We then focused in the problem of finding non-trivial (hairy) AFSSS black hole solutions where the Ricci scalar is 
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FIG. 9: (color online). Examples of numerical solutions for the Starobinsky Model 5 with q = 2 and As = 1. Left panel: Ricci 
scalar for three values of Rh (including Rh = 0 leading to the trivial solution R(r) = 0). R{r) —>■ 0 asymptotically, however, for 
the non-trivial solutions (dotted lines) the mass function M{r) does not converge (see the middle panel). Middle panel: mass 
function M(r) associated with the solution of R with Rh = 0.1 as shown in the left panel. Similar plots for M{r) are found, 
but not depicted, for the other non-trivial solution of R, while M{r) = const = r(,/2 when R{r) = 0, corresponding to the 
Schwarzschild solution. Right panel: metric components associated with the solution of middle panel. The metric components 
—gtt and Qrr and their product —gtt x grr = oscillate as r —>■ -|-oo corroborating that the resulting spacetime is not AF. 

Here R is given in units of Rst, and M and r in units of 1 / \/Rst {Go = c = 1). 
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FIG. 10: (color online). Example of a numerical solution for the Starobinsky Model 5 with g = 4 and A^ = 1. Left panel: the 
numerical solution shows that the spacetime is not AE as R{r) does not vanish asymptotically. It has rather a spurious De Sitter 
behavior when R{r) —>■ const > 0 asymptotically. Middle panel: the mass function M (r) does not converge asymptotically but 
decreases to a very large negative values. This is even opposite to the growing behavior M(r) ~ r® > 0 that should be expected 
if the spacetime were genuinely asymptotically De Sitter. Right panel: the metric components —gtt, grr and their product 
-gtt X grr = vanish asymptotically. This behavior confirms that the spacetime is not even genuinely asymptotically De 

Sitter where a cosmological horizon is expected at > th where gtt{rh) = 0 and grr{rh) = -too. 


not constant in the domain of outer communication of the black hole but varies with the radial coordinate r. Within 
that aim we provided the equations to study this scenario for an arbitrary f{R) model and derived the conditions for 
a regular black hole. We then proceeded to analyze some specihc models. Prior to a thorough numerical analysis, we 
studied the models under the scalar-tensor approach within the Einstein-frame, stressing that in vacuum f{R) gravity 
takes the same form as in the Einstein-scalar-field system. Therefore, this method makes possible to check if the 
available no-hair theorems for AFSSS can be applied, and thus, sparing an unnecessary numerical effort. Thus, for 
the models where the NHT’s apply, we concluded that geometric hair is absent. In those situations the numerically 
analysis simply confirms the NHT’s. 

For the cases where the resulting scalar-field potential ^ {(j)) does not satisfy the condition {(j>) > 0 required by 
the NHT’s, we turned to a detailed numerical study. Our conclusion is that we did not hnd any such hair in any of 
the models where the NHT’s do not apply. 

We also discussed some exact solutions that seem to represent AF hairy black holes, however, we showed that in fact 
the r dependence in the Ricci scalar is due to the presence of a deficit angle. Therefore such solutions are not really 
AF. Similarly we report numerical solutions that are not genuinely AF as the the mass function never converges to the 
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FIG. 11: (color online). Numerical solution for the Hu-Sawicky Model 6 with n = 1 and and ci « 0.38 and C 2 = 2.01 x 10 ® 
The units are like in Figure [To| but with Rhs instead of Rst- The description of panels is similar to Figure]^ 


Komar mass, a feature necessary to prove asymptotic flatness. Some of the solutions have an asymptotic oscillatory 
behavior but without encountering any singularity, while others produce singularities in the equations at a hnite r. 

It remains thus an open question to determine if hair exist in f{R) gravity. In order to settle the question in the 
affirmative it is tantalizing to depart from the Einstein-t/) system with a potential that has negative branches and 
that allow for hairy black holes 13 [IQ], and then perform an “inverse” conformal transformation to obtain an f{R) 
model. However, in practice this seems to be difficult in a closed explicit form. Moreover, even if this is possible, the 
resulting f{R) model would be rather artificial and would need to be submitted to the usual test (i.e. cosmological, 
Solar System, binary pulsar, etc.) to be better motivated physically, irrespective of the issue of hair. Some of the 
exact black hole solutions that have been reported in the literature (including in general relativity with exotic energy- 
momentum tensors) have been obtained using similar “tricks” or ad hoc confections, but presented afterwards in the 
more logical direction within the aim of enhance their merit. This is not the exception in f{R) gravity, where an 
f{R) model can be deduced by demanding the existence of some kind of exact solution, most of time, deprived of 
actual physical interest. For instance, one can impose grr = —^/gtt (in area coordinates), which is valid provided 
Tgfj j and then obtain a differential equation for f{R) which can be solved if the assumed form for < 7*4 (r) is 

simple enough. In fact, by using this method we were able to recover the model f{R) = f{R) = 2ay/R— a given by 
(19). Notice however that such solution was not AF, and thus, it cannot be used as a counterexample to the no-hair 
conjecture in f{R) gravity. 

In this paper we limited ourselves to the case of AF spacetimes. In a future investigation we will analyze the case of 
hair in asymptotically De Sitter and Anti-De Sitter spacetimes. For such spacetimes the no-hair theorems for the AF 
scenario require amendments due the different asymptotic conditions, and the numerical treatment, although similar 
to the one presented in this work, is sufficiently different to require a detailed and separate analysis. 
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Appendix A: Conservation of the energy-momentum tensor of matter in f{R) gravity 


The simplest way to show the conservation equation V°'Tab = 0 in f{R) metric gravity is departing from Eq. ([^ 
written as: 


fnRab 2 ^f 9ab b^ a 9ab^') /r — ^^ab • 


(Al) 


One can use the commutation relation = ^b^afn since fa is a scalar function (i.e. in f{R) metric gravity 

we assume the absence of torsion: the l.h.s of Eq. (All must be a symmetric tensor). So, applying the covariant 
derivative V“ to Eq. we obtain 


- \fR9abV‘^R " (V^VfcVa - VbV’^Va) Ir = R^^^Tab , 

where □ = V“Vq and fn = df/dR were used. According to the Bianchi identities 

1 


’'(^Rab — 1^9abF^ = 0 


therefore the first and third terms of l.h.s of Eq. (A2) cancel out, yielding 


Furthermore 


- (V“VbVa - VfcV^Va) /fl = . 

(V“VbVa - VfeV^Va) fR = (VeVfe - VfeVe) Va/fl = 9""Rcb/^dfR 

= Rbd^'^fR ■ 


This result implies that the l.h.s of Eq. (A4| vanishes identically, leading then to 

= 0 . 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 


We can summarize this result as a theorem: in f{R) metric gravity (under the same basic axioms assumed in GR 
concerning the manifold and the metric) the generalized tensor 


obeys a generalized Bianchi identity 


■— fRRab ^2^ b 9ab^'}fR 


= 0 . 


(A7) 

(AS) 


Appendix B: Supplementary regularity conditions 


In principle the regularity conditions at the horizon (34), (35) together with 6{rh) = 6h, R{rh) = Rh, M{rh) = rh/2, 
seem enough data to solve the system of ordinary differential Eqs. ([^, (10), and ( [II| ). Nevertheless, in view of the 
method that we devised to solve these equations numerically, care must be taken when evaluating at the horizon the 
right-hand-side of Eqs. (|^ and (11). Therefore at this place we also need to know the values of R'fj and 5'^ at the 
horizon when starting the numerical integration. If we do not impose these supplementary regularity conditions a 
mild numerical error is made at rH^ but we want to avoid this and impose also the right conditions on R'fj and 5'^. 
As we mentioned in Sec. IV in order to find these two conditions we require first to obtain an expression for R"' and 
5”, then as before, develop the quantities around r//, and finally impose that R'" and 5" are finite at the horizon. We 
simply provide the final outcome of this process, and the regularity conditions turn out to be as follows: 


5 ' 


= 4r^< r^ffuR 


/(lO/ 


- 13i?/ii 



f(2/- 


RJr 


+ r^ffR 


10/-i? 


(l3/fl 




+ 2/3 [e/ + RfR(2r^R - 3 )] U /i^i^[ 2/^(3 - + r^] [( - 4IrV + 2fR(l3r^R - GS^y^RfR 


+ 2r2(7rV2 + i8/^|^^+3/)) 


(Bl) 
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(B2) 


Given that these two regularity conditions are quite involved, we performed minimal tests to prove its validity. Like 
in Sec. IV we used the exact solutions provided in Section III First we used the solution Eqs. (14)-(I7|, notably for 
the model f{R) = kR?, and then we considered the model (181 with the solution provided by Eqs. (20)-(22). In both 


solutions (5(r) = 0. This means that the r.h.s of Eq. (BI) must vanish in the two cases. We verified that indeed this 
happens. 

On the other hand, the exact solution R(r) = Ri = const obtained from the model f{R) = kR ^ leads simply to 
= 0, whereas = 6/r^ for the model (18). In both cases we checked that the r.h.s of Eq. (B2) gives respectively 
these two values at the horizon. We are therefore confident that our expressions are correct and so we enforced them 


in the numerical treatment presented in Section V B 


Appendix C: Properties of the STT approach to f{R) gravity 


When introducing the action (36) in this way, we can recognize two things: 1) let us consider the functions 


H{R,Q) = f'{Q){R-Q) + f{Q) = Rf{Q)-[Qf'{Q)-f{Q)], and H{R,Q,x) = Rf'{Q)-N{Q,x), where N{Q,x) := 

Qx — fiQ), and in this Appendix a prime indicates differentiation with respect to the argument of the corresponding 
function. So, assuming f{Q) to be a convex function, i.e., f"{Q) > 0, then clearly BqN = 0 if x = f'{Q), and the 
inverse of this function Q(x) allows to define L(x) := N{Q{x)^x) = Q{x)x ~ fiQix))- la ttiis way we see that L{x) 
is no other than the Legendre transformation of f{Q), and H{R,Q{x)) = H{R,Q{x),x) = Rx ~ L{x)- Moreover, 
H{R,Q{x)) defines in turn the Legendre transformation of L{x) (at least in the region where L''{x) = Q'(x) > 0) 
which in this case the condition d^H = R — L'(x) = R — Q{x) = 0; simply leads to i? = Q{x)- 2) The condition for 
the second Legendre transformation can be imposed in the action by considering x = f'iQ) as a Lagrange multiplier, 
so that the variation with respect to x leads to R = Q. This is the formal construction when treating f{R) theories 
as STT, and in practice it is achieved by taking the action (37). 


In the following we perform explicitly the transformation between the original f{R) theory and the special class of 
Brans-Dicke model wbd = 0 supplemented with a potential. 

The Brans-Dicke action with a potential W (4)) is given bv|57j 


lBD[gab, ^ 


f 






t [dabj '0] 


(Cl) 


For comparison with the f{R) model we shall focus only on the case wbd = 0. Thus, the held equations read [57] 

W 

^Gab — biRab T (V4^ ^^9ab ^ 

□ 4 > = ^(^KT + <l>d^W -2w'j . 

On the other hand when introducing x = Ir Eq. ([^ simply reads 

XRab - ^Jgab - (VaVb - Pab^) X = f^^ab (C4) 


(C2) 

(C3) 


where /(x) = f{R{x))- Moreover, this equation can be written as 


XGab 


uTab + (ya^b - ffa&n) X “ ^9ab X-R(x) “ f (Rix)) 


(C5) 


where we have made explicit the functional dependence R{x)i which means that if f"{R) > 0, one can in principle 
invert the dehnition x := fR{R), and obtain R{x), and thus /(x) = /(i?(x))- Therefore, if we choose the potential 
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W (<i>) = — /(i?($)) = U (<&) and identify $ = x, then Eq. (C2) becomes exactly Eq. (C5). Moreover, the trace 

of Eq. (|C5[) reads 


^ + 2 / - X-R) 


With the above identification of the fields $ = x ^-nd the potential W we appreciate that Eq. 


(C6) 


also becomes 


Eq. (C6| where one can easily verify that the expression 2/ — x^R coincides exactly with ^d,s>W — 2W. Henceforth, 
we cone ude that f{R) theory is equivalent to a Brans-Dicke theory with wbd = 0 and a potential U{^). 


Appendix D: Examples of AF spacetimes with a deficit angle 


SSS spacetimes with zero charge that are asymptotically flat except for a deficit angle have the asymptotic form 


ds'‘ 


(. A 2 M\ , 


dr^ 


+ r'^dn‘^ . 


1_ A_ 


(Dl) 


After a redefinition of coordinates r = (1 — t = (1 — A) and M = Ma_]jma{^ — A)^/^ the metric acquires 

the standard angle-deficit form 


1 - 


2Madma \ 




^ 2 ^ _ ^MaoA£a'J 


+ {i-Aydn^ . 


(D2) 


Under this parametrization the coefficient Madma is then identified with the ADM mass associated with this kind 
of spacetimes [22] • In the same way, an SSS metric which is ADS or AADS with a deficit angle 


ds‘^ 


( ^ 2M Ar^ \ 


dr 


dr^ 




_ A — — — —^ 


can be transformed into 


1 - 


2Madma Af ^ 


3 / 


dt 




dr 


_ '^Madma. _ Af ^ \ 


+ (l-A)PdD" . 


(D3) 


(D4) 


Notice that the cosmological constant A did not require to be redefined in order to obtain the standard metric (D41. 
We then need Madma = 0, A = 1/2 and A = A^o in order to recover the metric (19) in the standard form when 
Q = 0. 

Finally, the metric of an SSS that is ADS or AADS with a deficit angle and endowed with a charge q 


ds'^ 


, 2M Ar^ , 2 


dr^ 


- A- 


2M _ AH I H 
r 3 




is transformed into 


ds^-- 1 - 


2Madma Ar^ 



+ (1 - A)r‘^dn^ 


(D5) 


(D6) 


taking q = Q{1 — A). The quantity Q is presumably the actual charge when A ^ 0. Again, taking Madma = 0, 
A = 1/2 and A = Aqo we recover the metric (19) written in the standard form but now with Q ^ 0 given by 
Q = ±g2 = ±Q2/4. 
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